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We investigate entanglement properties of multipartite states under the influence of decoherence. 
We show that the lifetime of (distillable) entanglement for GHZ-type superposition states decreases 
with the size of the system, while for a class of other states -namely all graph states with constant 
degree- the lifetime is independent of the system size. We show that these results are largely 
independent of the specific decoherence model and are in particular valid for all models which deal 
with individual couplings of particles to independent environments, described by some quantum 
optical master equation of Lindblad form. For GHZ states, we derive analytic expressions for the 
lifetime of distillable entanglement and determine when the state becomes fully separable. For 
all graph states, we derive lower and upper bounds on the lifetime of entanglement. The lower 
bound is based on a specific distillation protocol, while upper bounds are obtained by showing that 
states resulting from decoherence in general become non-distillable or even separable after a finite 
time. This is done using different methods, namely (i) the map describing the decoherence process 
(e.g. the action of a thermal bath on the system) becomes entanglement breaking, (ii) the resulting 
state becomes separable and (iii) the partial transposition with respect to certain partitions becomes 
positive. To this aim, we establish a method to calculate the spectrum of the partial transposition for 
all mixed states which are diagonal in a graph-state basis. We also consider entanglement between 
different groups of particles and determine the corresponding lifetimes as well as the change of the 
kind of entanglement with time. This enables us to investigate the behavior of entanglement under 
re-scaling and in the limit of large number of particles N ^ oo. Finally we investigate the lifetime 
of encoded quantum superposition states and show that one can define an effective time in the 
encoded system which can be orders of magnitude smaller than the physical time. This provides an 
alternative view on quantum error correction and examples of states whose lifetime of entanglement 
(between groups of particles) in fact increases with the size of the system. 

PACS numbers: 03.67.-a, 03.65.Ud, 03.67.Mn, 03.65.Yz 



I. INTRODUCTION 



Since the early days of quantum mechanics, entangle- 
ment has remained at the focus of interest. While entan- 
gled states of microscopic samples of matter -such as a 
few atoms or ions in a trap- can nowadays be prepared 
and studied in the laboratory 0, the question whether 
entanglement can persist on a macroscopic (i.e. classi- 
cal) scale is still subject of an ongoing debate. The puz- 
zling consequences of such macroscopic entanglement - 
highlighted in the notorious gedanken experiment known 
as "Schrodinger's cat" 13] by Schrodinger in 1935- and 
the (as it seems) absence of entanglement in macroscopic 
objects and hence in our "classical world" suggests a 
mechanism which prevents the persistence of entangle- 
ment on a macroscopic scale. It is often argued that 
decoherence -i.e. interaction of a system with uncontrol- 
lable degrees of freedom of some environment- Q pro- 
vides a way to understand the absence of entanglement 
on a macroscopic scale. In particular, the decoherence 
rate is believed to grow linearly with the size of the sys- 
tem which would predict a rapid decay of entanglement 
for systems consisting of many particles. 

Such an argument can easily seen to be valid for certain 



entangled states, e.g. superposition states of the form 

|Gi/Z) = 1/V2(|0)®^ + |1)®^), (1) 

also called Greenberger-Horne-Zeilinger (GHZ) states 
which are states of N spins or qubits, that interact with 
uncontrollable degrees of freedom of the environment, 
e.g. described by a heatbath. GHZ states can be viewed 
as simple models of Schrodinger cat states and are in 
fact sometimes called cat-states. For GHZ states, one 
can indeed show that if k is the decoherence rate of a 
single qubit, then the rate at which the iV-qubit state 
decoheres is given by kN. However, the observation that 
multipartite entanglement becomes more fragile with the 
size of the system is valid for this specific state only, and 
a priori it is not clear whether a similar conclusion can 
be drawn for other multipartite entangled states. 

Moreover, the decoherence rate does not provide com- 
plete information about entanglement properties of a sys- 
tem. In the last few years a theory of entanglement has 
emerged, which allows for a more sophisticated and de- 
tailed investigation of the effect of decoherence on the en- 
tanglement properties of a multiparticle entangled state. 
Quite recently, we have shown in Ref. 4] that for GHZ 
states not only the decoherence rate shows a scaling be- 
havior with the system size, but also the lifetime of distil- 
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lable entanglement -that is the time after which (distil- 
lable) entanglement disappears from a system subjected 
to decoherence- in fact decreases with the number of 
particles N, confirming the previous reasoning. On the 
other hand, we have also shown in Ref. that for a 
class of other genuine multiparticle entangled states - 
most notable cluster states Q-, the lifetime of distillable 
entanglement does not depend on the number of parti- 
cles and thus the size of the system. This is in sharp 
contrast to the behavior of GHZ states and shows that 
genuine macroscopic entanglement can indeed persist on 
timescales which are independent of the size of the sys- 
tem. 

While the investigation in Ref. was limited to a spe- 
cific decoherence model corresponding in physical terms 
to individual coupling of particles to a thermal reservoir 
in the infinite temperature limit, we will show in this ar- 
ticle that these observations are largely independent of 
the specific model of decoherence. In particular, a simi- 
lar scaling behavior of the lifetime of iV-party distillable 
entanglement with the size of the system is obtained for 
all decoherence models dealing with a coupling of each 
particle to its own environment (or heat bath), e.g. de- 
scribed by a quantum optical master equation of Lind- 
blad form. The results can even be extended to collective 
(finite range) couplings of particles to the environment. 
While for GHZ states we provide analytic results for the 
lifetime of (distillable) entanglement, we calculate upper 
and lower bounds on the lifetime for all states which be- 
long to the family of graph states The lower bound 
is based on an explicit entanglement distillation protocol, 
while upper bounds are obtained by three different meth- 
ods. Using the first method we show that the completely 
positive map describing the decoherence process becomes 
entanglement breaking after a finite time. This im- 
plies that all initially entangled states become separable 
and thus the lifetime of all kinds of entanglement is finite. 
The second method is more specific to graph states and 
shows that graph states suffering from decoherence be- 
come separable after a finite time. This is done by using 
a dynamical description of graph states and by showing 
that the generating operations become separable. The 
third method is based on the partial transposition cri- 
terion and evaluates when the partial transposition with 
respect to a certain partition becomes positive. To this 
aim, we develop a method to calculate the spectrum of 
the partial transposed operator p^^ for any subset S of 
parties and all density operators p which are diagonal in 
a basis constituted by orthogonal graph states. 

We also consider entanglement between M groups of 
particles, i.e. partitions of the system into M parts. Each 
of the groups may consist of several particles, which are 
then considered as a single subsystem with a higher di- 
mensional state space. We analytically determine the 
lifetime of distillable entanglement between M groups of 
particles for arbitrary partitionings for GHZ states and 
again derive lower and upper bounds for all graph states. 
In this way we study the change of the kind of entangle- 



ment with time and e.g. show for GHZ states that the 
effective size of entanglement, i.e. the maximum number 
of entangled subsystems, decreases with time and entan- 
glement eventually becomes bipartite before it vanishes 
completely. If we associate a specific spatial distribution 
with the particles, e.g. spins distributed on a lattice, 
one can choose certain partitionings that correspond to 
a re-scaling of the size of the subsystem, as it is used 
in statical physics. We study in particular the behavior 
of distillable entanglement under coarsening of the parti- 
tions, that is under re-scaling of the size of the subsystem 
in the asymptotic limit N ~f oo. For cluster states (and 
all other graph states with constant degree) we show that 
the lifetime of distillable entanglement is largely inde- 
pendent of N and thus the same on all scales. For GHZ 
states, however, we find that whenever the size of the 
subsystems is finite, distillable entanglement vanishes af- 
ter an arbitrary short time on all scales. Only if the size 
of the subgroups become macroscopic themselves (in the 
sense that A'' systems are divided into a fixed number M 
of cells whose size N/M grows to infinity as ^ oo) the 
lifetime of distillable entanglement (between the M cells 
of macroscopic size) becomes finite and scales to leading 
order as 1/{kM). We also consider the lifetime of en- 
coded entangled states. When considering entanglement 
properties between groups of particles where each group 
constitutes a logical qubit of a (concatenated) quantum 
error correction code, one can define an effective time 
for the encoded system, that incorporates the error cor- 
rection procedure. The effective time can be orders of 
magnitude smaller than the physical time. In this way 
one can show that the lifetime of entanglement between 
groups of particles can even increase with the size of the 
system. 

The paper is organized as follows. In Sec.|n]we intro- 
duce decoherence models -most notable individual cou- 
pling of a single particle to a reservoir described by a 
quantum optical master equation of Lindblad form as 
well as Pauli channels- which we deal with throughout 
the paper. In Sec. Illll we introduce basic concepts of en- 
tanglement theory. In particular, we review the concepts 
of separability and distillability in multiparticle systems 
as well as the the partial transposition criterion. We also 
define lifetime of entanglement with respect to certain 
partitionings of the system. In Sec. IIVI we determine 
the lifetime of iV-party distillable entanglement of GHZ 
states for decoherence described by depolarizing chan- 
nels as well as general quantum optical channels. In 
Sec. we first review the concept of graph states 
in Sec. IV Al and then derive lower and upper bounds 
on the lifetime of A^-party distillable entanglement for 
graph states subjected to decoherence. We generalize our 
results to weighted graph states in Sec. IV Fl In Sec. IVII 
we consider entanglement between groups of particles for 
GHZ states and determine the lifetime of encoded entan- 
gled states in Sec. IVI CI Wc summarize and conclude 
in Sec. IVIII while some technical details e.g. regarding 
the partial transposition criterion and the correspond- 
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ing upper bound on the lifetime for mixed states which 
are diagonal in a graph-state basis can be found in the 
appendices. 



II. DECOHERENCE MODELS 

We consider a single two-level system (qubit) coupled 
to an environment which is described by a thermal reser- 
voir. The evolution of this qubit is governed by a general 
quantum optical master equation of Lindblad form 
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-i[H, p] + Cp, 



(2) 



where H describes the coherent evolution while incoher- 
ent processes are described by the superoperator C. We 
have 
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with a± = l/2{axiif7y) and 2C > B. While parameters 
B, C give the decay rate of inversion and polarization, 
s e [0, 1] depends on the temperature T of the bath. 
More precisely s = limt^oo ( )t, where s = 1/2 cor- 
responds to T = 00. It is straightforward to solve this 
master equation We consider the case H = 0, i.e. 
solely decoherence. The eigenoperators and correspond- 
ing eigenvalues of C can readily be determined and one 
finds 
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For p = p(0) = ^]l + a- (7 = (To+ dx^x + O-yf^y + [Oz 

(2s — l)/2](Tz, we have that 



p{t) = e^V(O) = o-Q + e ^"{axCTx + ayOy) 

+e-^*[az-(2s-l)/2]az. (8) 

Equivalcntly, one can describe the resulting completely 
positive map (CPM) Et with p[t) = £tp as follows: 
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£tP = ^Aj(i)crjpcrj- 

+p{t)[azpt + \pcfz — ioxpoy + idypox] (9) 

with 



Ao(t) = -(l + 2e-^* + e-^*) 



(10) 



Ai(t) 



A2(i)--(l-e 
1 



-Bt\ 



hit) = -(l-2e 



-ct 



4 

2s- 1 



(1-e-^*) 



(11) 
(12) 

(13) 



In Sec. IVCl we will discuss the entanglement properties 
of this map and show that it (except for some singular 
cases) becomes entanglement breaking after some finite 
time. For s = 1/2 and B = C = k, Eq.Q describes the 
coupling of the particle to a thermal bath in the large 
T limit equivalent to a so-called depolarizing channel 
(white noise): 



Vp^pit)p 
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.(14) 



For B = 0, C = K and arbitrary s, Eq.® describes 
instead the coupling of the particle to a reservoir, that is 
equivalent to a dephasing or phase flip channel: 



Vp^p{t)p 



l-p{t) 



(p + aspas) with p{t)^e '=*(15) 



Finally, choosing s = 1 and B = 2C = k, one obtains the 
decay channel (pure damping): 



Vp = EipEl + E2PEI 



(16) 



with the Kraus operators Ei 

V7' 
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. Here 7(t) = 1 



denotes the decay rate 



for the decay from level |1) into level |0). 
For a system consisting of several particles, we shall be 
interested in the effect of decoherence on the entangle- 
ment properties of this system. We consider as a deco- 
herence model individual coupling of each of the qubits 
to a thermal bath, where the evolution of the fc*^ qubit 
is described by the map £). given by Eq. with Pauli 
operators tJj acting on qubit k. We will be interested in 
the evolution of a given pure state \^) oi N qubits un- 
der this decoherence model. That is, the initial state j^*) 
suffers from decoherence and evolves in time to a mixed 
state p{t) given by 



p(t) = £l£2...fAr|*)(*| 



(17) 



Disregarding a physical description in terms of an un- 
derlying interaction between the system and its envi- 
ronment, we will in the following also consider decoher- 
ence due to individual noise processes of the particles 
described by some Pauli channel: 



Vp. 



3 3 
^Piit)o-ipa^ with C^Piit) = 1) 

i=0 i=0 



(18) 



These noise channels are of particular interest in quan- 
tum information theory, especially in the study of fault- 
tolerance of quantum computation. This class contains 
for example: 
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1. for po = and pi = P2 = P3 = the above 
depolarizing channel; 

2. for Pi ~ \i the quantum optical channel according 
to Eq. ® with ^ = (s = i); 

3. for pq ~ p^—p2 — ^ and = the above 
dephasing channel; 

4. for Po — P2 = P3 = and pi = the bitflip 
channel 

In the remainder of the paper, we will analyze the time 
dependence of the entanglement properties of the deco- 
hered state p{€) for different initial states \^). The depo- 
larizing channel is of particular interest, since the deco- 
hered state due to an arbitrary noise channel can be fur- 
ther depolarized to some state, which might also be ob- 
tained by some depolarizing channel. Moreover, among 
the stated noise models the depolarizing channel is the 
only channel, that is basis independent, i.e. invariant un- 
der unitary transformations. We will frequently use the 
Pauli channel and will describe the entanglement prop- 
erties of p(t) in terms of the parameters pi. Nevertheless 
one has to keep in mind, that the time dependence itself 
is already included in the parameters pi = Pi{t). 

III. SEPARABILITY, DISTILLABILITY AND 
LIFETIME OF TV-PARTY ENTANGLEMENT 

For the lifetime of entanglement it is not only necessary 
to specify the underlying decoherence model, but also the 
very notion of multi particle entanglement itself. This is 
mainly due to the fact that multi party entanglement is 
a subtle issue in quantum information theory (see e.g. 

Apart from some special cases, the existence of 
an entanglement measure, that is satisfying for informa- 
tion theoretic purposes as well as applicable and calcu- 
lable for mixed states, is still an open problem. In the 
following we will therefore concentrate on the discussion 
of two qualitative entanglement criteria. Throughout the 
paper we will consider N two-level systems (qubits) with 
corresponding Hilbert space Ti = ((T^)**^. The N par- 
ticles are distributed among N parties 1, . . . , A^. Start- 
ing with a pure GHZ or graph state we will consider in 
Sec. lIVI andlVlthe A^-party separability and distillability 
properties of the decohered state p{t) (see Eq. (|17|l '): 
On the one side of the scale the state p{t) can still be 
N -party distillahle entangled, as it is the case for the cor- 
responding pure states in question. Hereby we call p{t) 
iV-party distillable, if any other true A^-party entangled 
state 1$) can be obtained (distilled) asymptotically from 
multiple copies of p under local operations and classical 
communication (LOCC) [Tol.lTH: 

P^' -~>LOCC |$)($| . (19) 

We remark that in the multi-copy case all true A^-party 
entangled states are equivalent since they can be trans- 
formed into each other by LOCC. That is, the condition 



that any true A'' party entangled state can be created 
can be replaced by the condition that some A'^-party en- 
tangled state, e.g. the initial pure state, can be created. 
Disregarding the practicability of the underlying distil- 
lation protocol, the state p{t) is then as useful as any 
other entangled state and therefore can in principle be 
regarded as a universal resource for quantum informa- 
tion processing such as quantum communication. 

On the other end of the scale, p{t) might have also be- 
come completely separable or classical in the sense that it 
can be described by a classical mixture of product states, 
i.e. p is N -party separable, if 

p(t)=^PfepW®pf (20) 

k 

If a state is completely separable, it is no longer entangled 
whatsoever. In between these two extremal cases, p{t) 
can contain different types of blockwise entanglement, 
which we will discuss in more detail in Sec. EH There 
we will consider different partitionings of particles into 
M groups (M < A^), where each group forms a subsys- 
tem with a higher dimensional state space and consists of 
several particles. M-party distillability [separability] can 
then be defined with respect to a given partitioning in a 
similar way, where the notion of local operation has to be 
adapted accordingly. Moreover we will call p{t) M-party 
distillable, if there exists at least one partitioning, with 
respect to which p{t) is M-party distillable. 

Based on the notion of M-party separability and dis- 
tillability, one can define lifetime of entanglement. A 
A^-party state which is subjected to decoherence 

for time t evolves into a mixed state p(t). The lifetime 
of A^-party distillable entanglement is given by the time 
after which the state p{t) becomes non-TV-party distill- 
able. This implies that lower bounds on the lifetime of 
distillable entanglement can be obtained by showing that 
the state p{t) is distillable, while an upper bound can 
be obtained by proving non-distillability of p(t). When 
considering partitions of the system into M groups, the 
lifetime of M-party entanglement with respect to a given 
partition is defined accordingly. We refer to the lifetime 
of M-party entanglement as the time after which p{t) 
is non-distillable with respect to all M-party partitions. 
In a similar way, one can define a lifetime with respect 
to the separability properties of p{t). 

In order to determine entanglement properties of the 
mixed states in question, we will continuously make use 
of the partial transposition criterion jl2l IT^ , an entan- 
glement criterion which provides necessary conditions for 
distillability and separability. The partial transposition 
is defined for bipartite systems only, while a system can 
in general consist of several parties. Making use of the 
concept of partitionings of the system, in particular con- 
sidering all bipartitionings, one can use the partial trans- 
position criteria also for multipartite states. Let A de- 
note a subset of m parties ki, . . . ,km- In general, given 
an operator X acting onQ!'^^ we define the partial 

transpose of X with respect to the first subsystem in the 
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basis {|1), |2), . . . , Id^)}, X^-*, as follows: 

X^-EE ^ Y.{zMX\3,l) \:hk){iA- (21) 

i,j=l k,l=l 

A hermitian operator X has a non-positive [positive] par- 
tial transpose (NPT) [(PPT)] if X'^^ is not positive [pos- 
itive] respectively. That is, X'^^ is NPT if there exist 
some l^-) such that (^'IX^-*!^') < 0. 

The positivity of the operator p^^ gives a necessary cri- 
terion for separability, whereas the non-positivity of p^^ 
is necessary for the distillability of the density operator p. 
In particular, if a bipartite density operator is PPT, then 
it is certainly not distillable This implies ^3] that 
if a multiparticle density operator p is PPT with respect 
to at least one bipartite partition, then p is certainly 
not A^-party distillable. On the other hand, positivity 
of all bipartite partitions is a necessary condition for iV- 
party separability. In the case of two dimensional systems 
(g) the PPT [NPT] criterion is necessary and suf- 
ficient for separability [distillability] 13, 14]. A detailed 
discussion of the application of the partial transposition 
criteria to multipartite systems can be found in Ref . |llj] . 



IV. LIFETIME OF N-PARTY ENTANGLEMENT 
IN GHZ STATES 

We start by considering the lifetime of A-qubit GHZ 
states 

lGi7Z) = 1/V2(10)^^ + 11)®^). (22) 

These states are special examples of states, that maxi- 
mally violate multi-partite Bell inequalities 15]. GHZ 
states have also become an interesting resource for multi- 
party quantum communication, e.g. in the context of 
secret sharing and secure function evaluation More- 
over they can be used to improve frequency standards 
[l7| . For the class of A^-party GHZ states the lifetime of 
most of the above entanglement properties can be deter- 
mined analytically. 



A. Large T limit of reservoir 

We start by reviewing the results of and consider a 
model of decoherence with individual coupling of each of 
the particles to a thermal reservoir in the large T limit. 
The process for a single qubit is described by Eq. (|14|l and 
corresponds to white noise with time dependent parame- 
ter p = p{t) — e^*** where k is a coupling constant. That 
is, we consider the state p{t) given by Eq. ((T7jl where 
the CPM £k is given by the the depolarizing channel Vk 
(Eq. It is straightforward to evaluate the effect 

of decoherence on this kind of states [Til ITsI ] . One finds 



that \GHZ) evolves to a state p{t) 
1 

p{t) = ^ Afc,...fc„Pfc,...fc„+M4'^ + ^^^), (23) 

with 

Pkik2...kM = \k1k2 . ■ ■kN){kik2 . . .k^l, (24) 

and cr± = (cr^ ±icry)/2, i.e. ct|^ = [00 . . . 0)(11 . . . 1]. 
It turns out that the coefficients Xkik2...kN fulfill 
Afcifc^.-.fcjv = •^fcifc2...fciv' w^ere kj = I - kj. In addi- 
tion, \kik2...kN depend only on fc = X^jLi^j' that is 

Afcifc2...fc« = h,i2...iM = Afc if kj = Y,^^ I J = k. This 
implies that = XN-k and one finds 

N 

^ - —■ (25) 

States of the form Eq. with Xk^.-.k^ = Xy.^ j.^ can 
equivalently be written as density operators which are 
diagonal in a basis consisting of orthogonal GHZ states. 
Such density operators have been completely character- 
ized with respect to their entanglement properties in 
Ref. In particular, it was shown that these states 

are A^-party distillable [separable] if and only if the par- 
tial transpose with respect to all possible partitions is 
non-positive [positive], respectively. One readily finds 
that the partial transposition with respect to a group 
which contains exactly k parties is positive, p(t)^^k > 0, 
if and only if ^3] \j? < XkXN-k, i-e. 

< 2Afe. (26) 

Making use of the fact that 

Ai > A2 > ... > A[^v/2], (27) 

it is now straightforward to determine the lifetime of dis- 
tillable A^-party entanglement as well as the time when 
the state becomes fully separable. From Eqs. H261 12 7|) 
follows that the lifetime of distillable A^-party entangle- 
ment is determined by Eq. H26I) with k = 1, as the par- 
tial transpose with respect to the partition one party - 
(A^ — 1) parties is the first one to become positive. Simi- 
larly, Eq. (|26|l with k = [N/2] determines the time after 
which the state p{t) becomes fully separable, as the par- 
tial transposition with respect to the partition N/2 — N/2 
parties is the last one to become positive. 

One observes that the critical value pcrit = e"***"^"', 
at which the partial transposition with respect to one 
party becomes positive increases with A^. This implies 
that for t > fcrit = the state is no longer A^-party 
distillable entangled and thus the lifetime r of true A^- 
party entanglement decreases with the size of the system 
as expected (see Fig. ^ and 15} . Note that finding the 
threshold value pcrit for a given A^ exactly is equivalent to 
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finding the roots of a polynomial of degree N (which can 
be done efficiently numerically) . One can obtain analytic 
upper and lower bounds on pcrit by approximating Afc by 
(l+p)W-fc(l_p)fc]/2A'+i or 2(l+p)^-'=(l-_£)^]/2^+i 
respectively, which is done explicitly in Sec. IVI Bl 
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FIG. 1: Upper bound on lifetime kt of A'^-party entangle- 
ment. 
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FIG. 2: Same as Fig. 1 but with double-logarithmic axis. 
Note that the same figures also display an upper bound on the 
lifetime kt of M-party entanglement in systems with N ^ oo 
particles for different M, as discussed in Sec. IVII In this 
case the numbers on the x-axis have to be considered as the 
number M of M-party entanglement in question. 



B. Arbitrary individual coupling to the 
environment 

We will now investigate the lifetime of distillable en- 
tanglement for more general decoherence models. While 
we continue to assume an individual coupling of parti- 
cles to independent environments -an assumption which 



is particularly well fulfilled if the entangled states in ques- 
tion are distributed among several parties-, we consider 
now couplings which are described by arbitrary quantum 
optical master equations of Lindblad form (see Sec. ITl|l . 
These models include as particular instances decay chan- 
nels, phase fiip channels and depolarizing channels. We 
consider the influence of decoherence -described by the 
CPM Eq. ®- on the GHZ state of N particles, i.e. 
the entanglement properties of the density operator p(t) 
which is given by Eq. IjlTI) . We use that one can write 
p{0) = \GHZ){GHZ\ as 



1 



P(0)= 3(^0 



(28) 



where Pq 



1-1- 



(cTj, ± icr j,)/2. It is not difficult to see that the action of 
the map £ (Eq. ©) is given by 



£Po = aPo + {l-a)Pu 
£Pi = (l-c)Po + cPi, 
£a± = ba±, 



(29) 



where we introduced the new variables a, 6, c which are 
given by 



a = s + (1 — s)e 



-ct 



(1 - s) + se 



-Bt 



-Bt 



(30) 



It is now straightforward to determine the action of the 
map £±£2 ■ ■ ■£n on the state p{0). One finds that the 
resulting density operator is of the form Eq. 123(1 . The 



coefficients A 
Afc 



only depend on ^ 



N 



^ kj , where 



1 
2 

2 



(c'-^l - c)^-'^ + (1 - a)'^a^-'=) 



(31) 



The condition that the partial transposition with respect 
to k parties is positive, p'^^k > reads 



K/2V < AfcA 



(32) 



We remark that in contrast to the discussion in Sec. llV Al 
here we have Afc ^ ^N-k- This means that non-positive 
[positive] partial transposition with respect to all parti- 
tions is no longer a sufficient condition for A'^-party distil- 
lability [separability] respectively. However, one can still 
use the partial transposition criterion to obtain lower and 
upper bounds on the lifetime of distillable entanglement. 
In particular, if the partial transposition with respect to 
at least one partition is positive, then the state p{t) is 
certainly no longer A-party distillable. 

To obtain an upper bound on the lifetime of GHZ 
states, we make use of the following facts: (i) Xk^N-k > 
(ac)^-'=[(l - a)(l - c)]V4 and (ii) exp{-NBt)/4 > 
[6^/2]^. While (i) can be checked by direct computa- 
tion, (ii) follows from 2C — B > (see Sec. 1^). Using 
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(i) and (ii) together with Eq. H32|l . one obtains that p{t) 
certainly has positive partial transposition with respect 
to any group of k parties if 



N > k 



log(Qc)-log[(l-a)(l-c)] 
log (ac) + Bt 



(33) 



provided that s ^ 0, 1 and B > 0. We remark that the 
(singular) case s = corresponds to a decay channel, and 
for such a channel we have that the state p{t) has non- 
positive partial transposition for all times t. Whenever 
the temperature of the bath is however not zero (i.e. s ^ 
0,1) we have that for any time t there exists a finite 
number A'o (given by the right hand side of Eq. H33() 
with k = 1) such that for N > Nq particles the state p{t) 
is certainly no longer distillable. Thus we have -as in 
the case of depolarizing channels- a scaling of the (upper 
bound on) lifetime of distillable entanglement with the 
number of particles N. If is sufficiently large, the 
(upper bound) on the lifetime goes to zero. 

A lower bound on the lifetime of p{t) can be obtained 
as follows: We have that a state of the form Eq. H23() 
can be depolarized by means of a (stochastic) sequence 
of local operations and classical communication (see Ref. 
[llj l') such that the resulting state has new coefficients 
Afcifc2...fc„ which fulfiU 



Afeife2.--fc]V + 



kik2 



(34) 



and hence Afc = \N-k- It follows that the depolarized 
state p(t) is distillable if 



r/2] > Afc, 



(35) 



for all k. 



One can upper bound Afc by Afc < AJ, 

\N-k ^N-kt-i _^fc „fc/-| „\N-k „N-ki 



max(a'=(l - a)^-^ a^-'=(l - a)^ c'=(l - c)^-^ c^-''^(l 
c)'') and obtains that p{t) is distillable if > 2Afc. 
By taking the logarithm of this equation, one obtains a 
bound on the number of particles N such that the state 
remains distillable for a time t. 



V. LIFETIME OF N-PARTY ENTANGLEMENT 
IN GRAPH STATES 

In the previous section the class of generalized GHZ 
states was shown to have a lifetime of entanglement, that 
decreases (except in some singular cases) with the num- 
ber of particles N in the system. We will now discuss 
the lifetime of iV-party entanglement in graph states and 
show, that for a significant subclass such as the cluster 
states the lifetime of distillable entanglement is essen- 
tially independent of N. After recalling some basic defi- 
nitions and notations, we will first derive a lower bound 
to iV-party distillable entanglement by providing an ex- 
plicit distillation protocol. We will then use three dif- 
ferent techniques to establish upper bounds to the life- 
time of iV-party entanglement. These methods apply 



to different decoherence processes and are interesting in 
their own, since they might find applications also in other 
problems not directly related to lifetime of states under 
decoherence. Finally we will extend our results to a more 
general class of so called weighted graph states. 



A. Basic definitions and examples 

Graph states are multi-particle spin states of dis- 
tributed quantum systems with interesting applications 
in quantum information theory: Special instances of 
graph states are codewords of quantum error correct- 
ing codes, which protect quantum states against de- 
coherence in quantum computation. Up to local uni- 
taries all stabilizer states can be represented as graph 
states 19j. For example, the CSS-codes correspond to 
the class of so called 2-colorable graphs . For this class 
of graph states entanglement purification procedures are 
known Q . These protocols even work in the case of noisy 
local control operations. Finally the class of cluster states 
arc known to be a universal resource for quantum com- 
putation in the one-way quantum computer |2lj| . For the 
study of genuine multi-partite entanglement graph states 
are particularly useful, since they allow for an efficient 
description even in the regime of many parties: Thereby 
the graph essentially encodes an interaction pattern be- 
tween the particles. Let G = {V, E) be a graph, which is 
a set of N vertices k &V connected by edges {k, 1} £ E 
that specify the neighborhood relation between the ver- 
tices. Starting from the state |+)^'' := 0fcgy 1+)'-''' , 
where |-|-) = --i= (|0) -I- |1)) denotes the eigenstate of ax 

with eigenvalue -1-1, the graph state \G) is obtained by 
applying a sequence of Ising-type interactions 



U, 



kl 



,-»f (i('"-^<'=))»(i<"-<t(')) 



(36) 



according to the interaction pattern specified by the 
graph, i.e. 



\G) 



n 



(37) 



{k,l}eE 



Graph states occur e.g. as a result of the Ising interaction 
between neighboring spins on a lattice after a specific in- 
teraction time . An example for a realization of such a 
system is based on neutral atoms in optical lattices [2^ . 
Alternatively graph states can be specified in terms of 
their stabihzer: For this let Nk = {l€V\ {k, 1} G E} de- 
note the set of neighbors of k. Then the graph state \G) 
is the unique state in (C^)®^, that is the common eigen- 
state to the set of independent commuting observables: 



9 = ^(fc) 



r(0 



(38) 



where the eigenvalues to all fc G are 1. The stabilizer 
Sg of the state is thus generated by the set {K*^ \ k e V}, 
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which imphes 

\G){G\=J2a. (39) 

In order to obtain a complete basis for (C^)®^ we will 
also consider the eigenstates \U)g = f^lG) of Kjf ac- 
cording to different eigenvalues Uk, i-e. 

K^\U)g = (-1)^'= \U)g ■ (40) 

Here and in the following, sets U C T/ as an upper 
index for operators will label those vertices where the 
operator acts non-trivially, e.g. 

-^f = (g) '^i'^ ■ (41) 
keu 

Moreover we will denote sets U and their corresponding 
binary vectors U = {Uk)kGV — {Ui, ■ ■ ■ , Un) over F^(the 
integer field modulo 2) with the same symbol. Finally k 
will also denote both the vertex and the corresponding 
one-element set {fc}. In this notation the stabilizer gen- 
erators can be written as = (J^a^'' and the original 
graph state is just that with an error syndrome corre- 
sponding to the empty set 0, i.e. \G) = \0)g- This is also 
notationally advantageous, since we will use both set and 
binary operations: E.g. for A, B ^ ^(^) — '^^ '^il^ 
write AUB, AnB and A\ B (A = V\ A) for the union, 
intersection and difference (complement) as well as A+B 
and (A, B) for the addition and the scalar product mod- 
ulo 2. The neighborhood relation in a graph is also often 
represented in terms of its adjacency matrix F: 

~ \ otherwise. ^^'^> 

In the spirit of the above notation we can therefore also 

write = ai'^ Uievi'^i'^r^ = '^^-^ = ^.'^P- 

Coming to some examples, we first note, that the class 
of multi-party GHZ states in Sec. IIVI is contained in the 
class of graph states, since the GHZ state in Ea. H22|) can 
be transformed by local unitaries into graph states cor- 
responding to the graphs depicted in Fig. O When con- 




FIG. 3: The graph states corresponding to the complete 
graph and the star graph are equivalent to the GHZ state in 
Eq. i'2'21 up to some local unitaries Q. 

sidering decoherence of a locally equivalent state, we re- 
mark that the underlying noise process has to be adapted 



according to the local unitary transformation. From this 
point of view the depolarizing channel in Eq. (|14|1 has the 
advantage that it is invariant under local unitary trans- 
formation and hence is basis independent. In the follow- 
ing we will also consider the class of cluster states in 1, 2 
or 3 dimensions (see Fig. 0)), which are of particular in- 
terest in the context of 'one-way' quantum computation 
j2l| | . For more examples and a discussion of equivalence 
classes of graph states under local unitaries and/or graph 
isomorphics we refer to 0, |23 . We will now discuss the 




FIG. 4: Cluster states in 2 and 3 dimensions form a universal 
resource for quantum computation in the framework of the 
one-way quantum computer (gj- 

entanglement properties from Sec. IIIII for the state p{t) 
(see Eq. (|17|l '). that is produced by different decoherence 
processes described in Sec. ^acting on different types of 
graph states. 

B. Lower Bound: An explicit distillation protocol 

We establish a lower bound on the lifetime of graph 
states by considering an explicit distillation protocol. In 
order to show that a mixed state p{t) is A^-party distil- 
lable, it is sufficient to show that maximally entangled 
pairs shared between any pair of neighboring parties can 
be distilled. This is due to the fact that these pairs can 
be combined by means of local operations (e.g. by tele- 
portation) to create an arbitrary A'^-party entangled pure 
state. We emphasize that we use the distillation of neigh- 
boring pairs only as a tool to show A^~party distillabil- 
ity. This does, however, not imply that the entangle- 
ment contained in the cluster state was in some sense 
only "bipartite" . One could in principle also use direct 
multi-party entanglement purification protocols, e.g. the 
one introduced in Ref. [2^] , however the conditions under 
which these protocols are applicable are in general more 
complicated to determine. 

First we will consider the case of decoherence of the 
particles due to the same individual Pauli channel Vp = 
Y^^=oPi{^)'^iP'^i ^^'^ ^6 ^ill show then how to extend 
these results to more general decoherence models. We 
will essentially follow the ideas used in in which the 
corresponding result was shown for the case of a depo- 
larizing channel, and make use of the following facts: 

(i) Measuring all but two neighboring particles, say fc, I 
of a graph state \G) in the eigenbasis of cr^ results in the 
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creation of another graph state with only a single edge 
{k,l} jBI. That is, the resulting state of particles k,l 
is up to local <Tz operations equivalent to a maximally 
entangled state of the form 

|$)^-^(|0).|0). + |1).|1).) , (43) 

where \i)x [\i)z] denote eigenstates of [ctz] respectively. 

(ii) The action of a Pauli channel Vk acting on particle 
fc of a graph state can equivalently be described by a 
map Mk whose Kraus operators only contain products 
of Pauli matrices and the identity, where here may 
act on particle k and its neighbors, i.e. particles which 
are (in the corresponding graph) connected by edges to 
particle k. 

Observation (ii) follows from the fact that al\U)G = 
{—l)^'alK^\U)G, where criKj' is an operator which 
contains only products of <7z operators at neighboring 
particles of particle j, and the identity otherwise. Sim- 
ilarly, the action of cr^ on graph states is up to a phase 
factor equivalent to the action of an operator which con- 
tains only products of operators acting on particle j 
and all its neighbors. That is, 

Vk\u)G{u\ = j2pdt)'^'\u)G{uw^ 

3 

= M,\U)g{U\ = Y.M*)S'\U)o{U\S^, (44) 




FIG. 5: Ring with seven qubits. 

the case of a ring (Fig. |SJ, for example, we have Sq = 

1,S', = at'^af+'\s!l = and 5| = 

(fc) 

"3 ■ 

We now apply (i), (ii) to establish a sufficient condi- 
tion when bipartite entanglement between neighboring 
particles can be distilled from the state 

p[t) = MiM2...MN\G){G\. (45) 

This allows us to obtain a lower bound on the lifetime of 
graph states. We concentrate on two specific neighboring 
particles, say k and I. One performs measurements in the 
eigenbasis of Uz on all but particles k and I (We remark 
that measurements on all neighboring particles of parti- 
cles fc, I would also be sufhcient). It follows from (i) and 



(ii) that these measurements commute with the action 
of the CPM M1M2 ■■ - Mn on the graph state (which 
equivalently describes the action of Pauli channels on 
these states). That is, the resulting state after the mea- 
surements is given by M1M2 ■ ■ ■ MN\^)kj{^\ ® lx)(xl: 
where |x) is a state of the remaining {N — 2) particles, 
and 1$) is a maximally entangled state equivalent up to 
Gz operations (which can be determined from the spe- 
cific measurement outcomes) to |$) (see Eq. (|43ll '). We 
emphasize that the operator Mj only acts non-trivially 
on particle j and its neighbors. This is due to the fact 
(see (11)) that the operators S*/, I — 0,1,2,3 -and thus 
the map M.j~ only effect particle j and/or its neighbors. 
It follows that in order to determine the reduced density 
operator of two neighboring particles Pki{t), one 

has to consider only the action of maps Mj which act 
on particles fc,/ or neighbors of fc or Z on the maximally 
entangled state |$), i.e. 

Pki{t)^{\{Mm{^l (46) 

where / = Nk U Ni U k U I. We have that the reduced 
density operator pki (t) is distillable if and only if its par- 
tial transposition is non-positive i.e. Pki{t)^'' It 0- 
To obtain a lower bound on the time until which the A'^- 
particle state p(t) remains distillable one has to consider 
all neighboring pairs determine the correspond- 

ing threshold value on the lifetime of distillable entan- 
glement kI^^ and take the minimum over all neighboring 
pairs {fc,Z} G E [S^. For graphs corresponding to pe- 
riodic structures (e.g. some lattice geometry), such a 
minimization is however not required. 

Thus we have that the threshold value is a function 
of the local degree (i.e. the number of neighbors) of the 
graph, but is independent of the number of particles iV. 
Note, however, that the degree of the graph may itself 
depend on iV -as it is e.g. the case for GHZ states-, 
which then implies that the threshold value will indeed 
depend on N . In all cases where the degree of the graph 
is independent of N -which is e.g. the case for all graphs 
corresponding to some lattice geometry, such as 2D/3D 
cluster state, hexagonal lattices, lattices with finite range 
interactions, etc.-, we have no scaling with TV, i.e. the 
lower bound on the lifetime of entanglement is indepen- 
dent of the number of particles N . These results can also 
be understood in the following way: The measurement 
in the neighborhood of particles k and I disconnect these 
two particles from the remaining system, which implies 
that errors occurring in some outside area do not influ- 
ence the two particles in question. This insight is also 
used in the following Sections and allows one to show 
that the behavior of cluster states is not a consequence 
of the specific decoherence model but rather a general 
feature of such states. 

The exact dependence of the distillability properties of 
Pfei(i) (and thus the threshold value p^') on the graph G 
can be determined as follows: For j S \ Ni the action 
of A4j can be described by a phase-flip channel acting 



10 



solely on particle k, where a phase flip channel acting on 
particle k is defined by 

Mfp ^ p.p + + (47) 

and we find — 1 — 2(pi + P2) ~ 2(po + Pa) - 1- The 
action of A4j for j G Ni\ Nk can similarly be replaced 
by a phase-flip channel acting only on I. Moreover the 
action of A4j when particle j E Nk H A^; is a common 
neighbor of particles fc, Z is given by a correlated phase- 
flip channel, 

V = + ^-^ip + (48) 

where pzz = 1 ^ 2(pi +^2)- Note that the sequential 
application of each of these channels, say the correlated 
phase-flip channel with parameter for \NknNi \ times, 
is equivalent to a single application of the same chan- 
nel with new parameter p = pl^^nWil^ Finally the Pauli 
channels Aik and A4i have also to be taken into account. 
In any case the resulting state pki is diagonal in the "Bell- 
basis" { I , Ictz I $) , (Tz 1 1 $) , (Tz cr^ I $) } , where (p is given by 
Eq. H43() . One can now easily determine pki{t) for any 
graph G and thus the condition when pki{t) (Eq. (j^ l 
has non-positive partial transposition and is thus distil- 
lable. After some algebra, one obtains that the stated 
protocol yields distillable entanglement between k and I 
if 

• for the depolarizing channel Dp ~ pp + {1 — p) ^1: 

p\N,\ + l+p\N.+N,\^^\N,\ + l (49) 

holds , 

• for the bit-flip channel "Dp = pp+ {p + <yxP<^x)'- 

p\N.\ _^ p\N,+N,\ ^ ^ ^ (50) 

holds, 

• for the phase-flip channel Vp = PP + 
^ (p + azpaz): 

V2-l<p<l (51) 

holds, 

• for the quantum optical channel with /i = 0, i.e. 
Pi = P2 ^ P and p3 = g for < p, g < i; 

(1 - 2{p + q)) ((1 - 4p)l^'=l + (1 - 4p)l^'l + 

(1 - 4p)l^'=+^'l-2 (1 _ 2{p + <?)))> 1 (52) 

holds. 

A lower bound on the lifetime of distillable entangle- 
ment under decoherence due to one of the above Pauli 
channels can then be derived by solving the correspond- 
ing polynomial inequalities. From Ea. H51|) it follows for 



example that in the case of the phase flip channel the 
lower bound obtained by this distillation protocol is the 
same for all graph states. This can be understood by the 
fact that here only the two individual dephasing channels 
acting on k and I (and not those of their neighbors) are 
relevant for the decoherence of the bell state |$) between 
k and I. For the the bitflip and the depolarizing chan- 
nel the critical value p< for p, which is proportional to 
the fldelity with the original pure graph state, increases 
with \Nk\, \Ni \ or |A^fc-|-iV;|. Similarly, the critical values 
for {p,q) decrease with \Nk\, \Ni\ or jiV^ -I- A^/|, since p 
and q now represent the error probabilities instead of the 
fldelity with original pure graph state. In the following 
we will consider the condition H49() for the depolarizing 
channel with p = e^'^* in more detail: In the case of the 



-log(p„i,) 




GHZ bound 

20 40 60 80 100 N 



FIG. 6: Under individual coupling due to the same depolar- 
izing channel the lower bounds on nt to the lifetime of dis- 
tillable A''-party entanglement for the ID- , 2D-, 3D-cluster 
state remain constant for arbitrary system sizes . For the 
A^-party GHZ state the lower bound as well as the exact value 
for Kt according to Eq. (I26I I. until which GHZ state remains 
distillable entangled, strictly decreases and goes to zero as 
N —f 00. 

A^-party GHZ state both representations of Fig.|31yield to 
the same polynomial inequality 2p^ + p^ > 1, since the 
depolarizing channel is invariant under local unitaries. 
This polynomial inequality can be further estimated from 
above giving Kt > = o.emi (jepjctg^i pjg^ |^ 

the corresponding critical value for Kt is indeed always 
below the exact value given in Eq. H26|l and decreases 
with the number of particles N. For linear chains or rings 
i\Nk\ = \Nk+i \ = 2, \Nk + Nk+i\ = 4) one flnds a thresh- 
old value p< — 0.7167, which gives a lower bound on the 
lifetime nt^ = 0.3331. That is, for p > p<: {nt < K,t<^ 
the state p{t) is certainly iV-party distillable using this 
speciflc protocol. For cluster states corresponding to a 
regular 2D [3D] lattice we have that neighboring parti- 
cles corresponding to inner vertices, i.e. with |A'^fc| — 
\Ni\=A,\Nk + Ni\=%[\Nk\ = \Ni\=&,\Nu + Ni\ = 12], 
will give the most sensitive polynomial inequality p^ (p^ + 
2) > 1 [/(p5 + 2) > 1]. One hence finds p< = 0.8281 
(Kt< = 0.1886) [p< = 0.8765 {nt^ = 0.1318)] respec- 
tively. As it can be seen in Fig.^l whereas the lifetime of 
iV-party GHZ states decreases with the size of the sys- 
tem, cluster states do not show such a scaling behavior. 
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since the derived lower bounds for cluster states remain 
constant. 

In the following we derive a more handy expression for 
the critical values for p and nt: For fixed degrees \Nk\ and 
\Ni\, one finds that the strongest lower bound on the life- 
time (which is thus also valid for all other configurations 
of this kind) is obtained for \Nk O Ni\ = 0. This can be 
understood as follows. Assume that for some given graph 
one changes the graph such that the degree of two neigh- 
boring vertices k, I increases by 1, i.e. \Nk\ = l-^fel + 1 and 
I TV; I = I A'/ 1 -I- 1. The first possibility is that this increase is 
due to the addition of a single common neighbor of k and 
I, i.e. \Nk n iV;| = |iVfc n 7V;| -I- 1, which leads to the con- 
dition for distillability pl^'^+^'l + ^1^^1+2 + p|w,|+2 > i. 
In the second possibility the neighborhood of both par- 
ticles k and I is increased by two different particles, i.e. 
\Nk + Ni\ = \Nk + Ni \ + 2. In this case one obtains the 

> 1 for distiU- 



\Nk\+2 



condition _|_ p 

ability. Clearly, the second condition will give a larger 
value on p and thus provides a stronger bound on the 
lifetime of distillablc entanglement. Intuitively, this can 
be understood from the fact that adding a single joint 
neighbor corresponds to a single additional noise chan- 
nel with correlated phase noise, while adding two inde- 
pendent neighbors corresponds to two independent noise 
channels acting on particles k and I . The influence of two 
independent noise channels is larger than of a single (cor- 
related) noise channel. In order to derive a lower bound 
we may therefore evaluate the polynomial inequalities 
of the different neighboring particles fc, Z, as if the value 
\Nk~\-Ni\ was maximal (i.e. \Nk{^Ni \ — 0), since this will 
give stronger/larger critical value than the critical value, 
which would be the solution to the exact polynomial in- 
equality . Under this simplification and by using, that 

p\Nk\ + l j^p\Ni\ + l > 2p(l^'=l + l^'l+2)/2 = 2p(l^fc+^'l+2)/2^ 

one finds that for 



p > 2" 



Vfcl + |JVjl+2; 



(53) 



the reduced density operator pki is certainly distillable. 
This leads to the lower bound on lifetime 



21n(i) 



\Nk 



\Nr\ 



(54) 



Taking \Nk\ and |iV/| to be the maximal degrees of two 
neighboring vertices in the graph, this leads to a universal 
lower bound for all graph states under depolarizing noise. 

We remark that the observed behavior, i.e. that the 
lifetime of multiparticle entanglement for cluster (and 
similar) states is essentially independent of the size of 
the system, also holds for more general decoherence mod- 
els. This follows from the fact that — similar to (ii) — 
the action of any CPM acting on graph states describing 
an arbitrary decoherence process can be estimated by a 
CPM whose Kraus operators only contain products of 
operators and the identity and thus the measurement 
(i) still commutes with the CPM. To this aim, we apply 
after the application of the CPM a local depolarization 



procedure which maps arbitrary density operators to op- 
erators diagonal in the graph state basis without chang- 
ing the diagonal elements [2J| . When restricted on graph 
states, the resulting action of the initial CPM (given by 
'^ki^k,iOkpOi, where Ok,Oi arc products of Pauli op- 
erators)) can then be described by a CPM specified by 
X]fe o,k,kOkpOl, where all operators Ok can be expressed 
in terms of products of cr^ operators. Only operators Ok 
which act non-trivially on particles k,l or their neighbors 
in the graph affect the resulting maximally entangled pair 
after the measurement (i), leading again to a threshold 
value which is independent of the size of the system for 
all those decoherence models where the number of such 
operators Ok is independent of N. This is for instance 
the case if each Ok acts non-trivially on a finite, local- 
ized number of subsystems. Therefore, the fact that all 
graph states with finite maximal degree, such as cluster 
states, the lifetime of distillable iV-party entanglement 
will remain finite, holds in particular for all decoherence 
models based on an arbitrary individual coupling to the 
environment described by a quantum optical channel in 
Eq. 0. In the following we will determine upper bounds 
on the distillable entanglement. 



C. Upper bound I: Noise operation becomes 
entanglement breaking 

In our first approach, we determine an upper bound 
to the lifetime of iV-party entanglement by considering 
the capability of the decoherence process to disentangle 
any state disregarding its specific form. Hence the upper 
bounds derived in this way will apply not only to graph 
states but to an arbitrary state. In turn this method will 
be restricted to coupling of the particles to individual 
environments described by an arbitrary channel of the 
form 



PljCTipCTj 



(55) 



like the completely positive map (CPM) in Eq.®. 
We now make use of the Jamiolkowski isomorphism 
between CP maps and states [l^: Let | = 

^ (|0)'=|0)'''' + |1)''|1)''') denote the maximally entan- 
gled state on system k and a copy k' of the system k. 
Then to each CPM V acting on particle k there uniquely 
corresponds a state 



kk — 'T^k 



(1$+) 



kk' 



(56) 



on the combined system of k and k' . The main fact which 
we will use in the following is, that 

(i) the CPM T> is entanglement breaking j7|, i.e. V'^p'^^ 
is separable for any (possibly entangled) state p''^ on 
the composite system, consisting of particle k and some 
other particles hold by the parties A, if and only if the 
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corresponding state a!^ is separable (with respect to the 
particles k and k'). 

Hereby the 'only if implication directly follows from the 
very definition of the map to be completely disentangling, 
whereas the other direction can be seen as follows: Given 
the state crl/^ one obtains the corresponding CPM via the 
inverse isomorphism, i.e. 

Vp''" = 2^ X tikko (|$+)*'-^-''($+| p'^" ® 4'=') , (57) 

where p is an arbitrary state on another copy /cq of system 
k, the projection onto 1$+)'^"'^ is performed with respect 
to the joint system kQ,k and V is now thought to map 
system ko onto system k' instead of k onto k. Now, if 
CTp*^ = p\® P2 is separable, then H57(l factorizes into 

fcfeo^$+|^fco ^ pfc|$+)'=fco ^ ^fc' (-53) 

disregarding, whetherp*^" ~ p'^"^ is itself entangled with 
some other parties A or not. The resulting state on sys- 
tem k' which corresponds to Vp is even independent of 
the input state p and thus cannot be entangled with the 
parties A whatsoever. 

In order to derive an upper bound to the entangle- 
ment of states suffering from decoherence due to individ- 
ual coupling V'' of the particles to the environment, one 
can determine the critical value for pij in l|55() , for which 
the state 

3 3 

i,j—0 

becomes separable and hence the CPM I?*"' become en- 
tanglement breaking. In (|59|l we have used the notation 
= cri|$+), where (|$o), |^'2), I'^s)) form a com- 
plete 'Bell'-Basis. In the following we will restrict to the 
same individual coupling V'^ —Voi the particles to the 
environment, which then only requires to test the sepa- 
rability of one state axi- In the case of Pauli channels 
Vp = Yl^i=o Pi ^iP'^i tliis task becomes particularly easy, 
since the state axi is diagonal in the above 'Bell' basis. 
Moreover for such Bell diagonal states the separability 
criterion reduces to the necessary and sufficient condi- 
tion, that all diagonal entries pi are smaller than 1/2, 
i.e. 

1 , , 

max Pi < — . (60) 

1=0,1,2,3 2 

This can be easily evaluated for the examples given in 
Sec.m For the depolarizing channel Vp = pp+{l—p) ^1 
the state p{t) has certainly become iV-party separable, if 

P < P> = • (61) 

Note that this condition provides a universal upper 
bound for all states exposed to individual depolarizing 
channels. For the quantum optical channel with p = 



in Eq. ® one arrives at the condition 2Ai + A3 > ^■ 
In the case of a general quantum optical channel with 
/i 7^ or an arbitrary noise channel of the form l|55(l one 
can instead use the fact, that for any two dimensional sys- 
tems k and k' the PPT [NPT] criterion, i.e. the positivity 
[non-positivity] of the partial transpose p^'^As necessary 
and sufficient for separability [distillability] [i3,[lJ]- Thus 
the CPM V is entanglement breaking, if and only if 

{pfY''>0. (62) 

For the general quantum optical channel © this leads 
after some algebra to the condition 

A? - p' > (^Ai + A3 - . (63) 

In terms of the original parameters B, C and s of the 
quantum optical master equation with the superoperator 
defined in Eq. (|3J| this inequality reads 

s(l-s) [e^*(l-e-^*)]'>l. (64) 

It is worth remarking, that in the terminology of quan- 
tum optics (reservoir theory) both the equilibrium value 
s and the decay rates B, C enter in the inequality in this 
multiplicative form. For the example of a decay channel 
i.e. A = Ai = A2 = M, we have > (A -I- A3 — ^) , which 
cannot be satisfied. Therefore the decay channel cannot 
become entanglement breaking and the multi party GHZ 
states are an example for states, that remain entangled 
under decoherence due to this channel (see Sec. IIV B|l . 
For the bitflip channel Vp = pp + {p + (JxpcTx) and 
the dephasing channel Vp = pp + {p + (JzPO'z) the 
upper bound p = obtained from Eq. (|60|l becomes triv- 
ial. 



D. Upper bound II: Noisy Ising interaction 
becomes separable 

In our second approach, we determine an upper bound 
ty on the lifetime of distillable entanglement by show- 
ing that after a certain time, the state p{t) becomes fully 
separable and is hence no longer entangled whatsoever. 
To this aim, we consider the (dynamical) description of 
graph states in terms of Ising interactions acting on a 
specific separable state. We determine the separability 
properties of the operator p{t) by considering the corre- 
sponding interactions which generate the state and show 
that for a given noise level, these operations itself become 
separable and hence are not capable of creating entan- 
glement. Consequently, also the state p{t) is separable in 
this case. The main advantage of this approach is that 
one does not have to consider the iV-particle state p itself 
and determine when it is fully separable (a task which is 
generally very difficult, especially if N is large), but has 
to consider only two -particle operations and determine 
when these operations are separable. 
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We make use of the following properties: 
(i) The graph state \G) corresponding to a graph G 
can be written as (see Sec. IV A|l 5] 



iG)= n f^*n+> 



(65) 



{k,l}i£E 



where Uki 



))/2«(l(')-.(0)/2 ^^^^ ^ 



1/V2(|0) + |1)). 

(ii) We will only take a^- noise into account. We 
thus restrict the following analysis to decoherence mod- 
els due to the same individual noise channel £, that can 
be decomposed into some noise channel £' acting af- 
ter a dephasing channel Vp = PzP + ^—^[p + crzP<Jz\, 
i.e. £ = £' o D. A more detailed analysis of the cases, 
for which such a decomposition is possible, is postponed 
to the appendix A. For the depolarizing channel I?fc(p) 
(Eq. H14I) with noise parameter p = e""*), such a decom- 
position is possible choosing pz — and 



£ kP = ^ P 



1 



P 



(66) 



This can be checked by direct calculation. 

We now investigate the influence of noise on the entan- 
glement generating unitary operation Uu and determine 
when the resulting CPM becomes separable. Since Uki 
commutes with 'Dj{pz), it follows that p(t) can be written 
as p{t) = £'{p{t)), where p{t) is obtained from the origi- 
nal graph state by considering only phase noise described 
by Vk{pz), i-e 



p{t)^\{V,{pz) n 



fci' 



(67) 



{k,l}eE 



Since p{t) is obtained from p{t) by means of separable 
operations, it is sufficient to determine the condition 
when p{t) becomes separable. In principle, one could 
also consider this additional noise to obtain a stronger 
upper bound on the lifetime, however the analysis be- 
comes more involved in this case as one has to deal with 
correlated noise. In the following we will therefore con- 
sider only noise resulting from phase-flip errors described 
by 'Dk, i-e. the map 



Vki{pz,qz)p = 'Dk{Pz)'Di{qz)UMpUli , 

for two different dephasing parameters Pz and qz 
this notation the operator p{t) can be written as 



p{t)= n 

(k,l)eE 



)!+>(- 



(68) 
With 



(69) 



For the vertex k with degree |iVfc| we have split up the 
action of the map Vkipz) into \Nk\ parts (one for each 
term in the product which involves Uki and thus particle 
k) by using a decomposition of the the map Vkipz) into 



T>k{Pz)p ■ 



U^kip'J" 

i=i 



(70) 



This leads to the parameter pV^'^"^ in^. If in Eq. ^ 
all maps 'Dki{pV^^''^ ,pV^^'^) at a fixed vertex k are sepa- 
rable, it immediately follows that also p{t) is fc-versus-rest 
separable since the following maps are local and act on a 
fc-versus-rest separable state. 

To determine the entanglement properties of 
'DkiiPzTqz), we make again use of the Jamiolkowski 
isomorphism between CPM and mixed states 

psf . In particular, we use that a CPM T) is separable 
and hence not able to generate entanglement if the 
corresponding mixed state D is separable, where 



Dk^k^Ul, = Vk,i.A'^+)k,kA^+\ ® I 1, (71) 



and separability has to be determined between parties 
(^1^2) and {lil2)- It turns out to be useful to define 



^(|6)feifej6),,,, + |i)fc,feji)/,zj. 



\4>ii) 



-^(|0)fcifejl);ii2 

^(|6)fcifej6),,,, 

■^(|0)felfe2|l)ili2 



|l)fclfc2|0)/ii2): 



|l)fclfc2|l)/li2): 



|i>fcife2|6)/,i2), (72) 



with 



/felfc2 
)felfc2 
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—{\m)i,i, + \n)i,i,), 

-i=(|00)z,,, -|ll),,zj. 



(73) 



One finds that the state Dk^k2iii2iPz,qz) corresponding 
to the map VkiiPz^qz) is given by 

1 

Dk,k2hhiPz,qz) = J2 (74) 

where Aoo = (1 + Pz){l + qz)/4:, Xoi = (1 + Pz){l - 
g.)/4,Aio - (l + (7.)(l-P.)/4,An = (1 - p,)(l - g,)/4. 
This state is separable with respect to {kik2) — (hh) if 
and only if Aoo ^ 1/2, as its partial transposition is pos- 
itive in this case. Note that for systems inO^^ posi- 
tivity of the partial transposition is a sufRcient condition 
for separability |13i . .14] . Although the system that we 
consider consists of two four-level systems, the resulting 
state has support only in a four dimensional subspace and 
thus the results about qubit systems can be directly ap- 
plied. We then obtain that the operator Dk^k2iii2{Pzi qz) 
is separable -and hence the CPM Vkiipz, qz) is separable 
and not capable to create entanglement- if and only if 



il+Pz){l + qz) <2. 



(75) 
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We now use the above result to obtain an upper bound 
for the lifetime of graph states under a decoherencc 
model, that obeys (ii). Due to Eq. (|75|l we find that 

the map 'Dki{pV^^''\pl^^^'^) is separable if 



(i+pyi^'=i)(i+pyi-'i)<2 



l/|JV,h 



(76) 



The threshold value p> such that state p{t) is fully sepa- 
rable is then obtained by considering all pairs of particles 
{k, I}, calculate the corresponding value and take the 
minimum over all {k,l}, i.e. p> = minpy. This ensures 
that all involved operators are separable for pz < p>- By 
estimating \Nk \ and |iV;| from above with maximal degree 
m in the graph, we arrive at the weaker upper bound 



P:. 



< 



{V2-iy 



(77) 



For the various decoherence processes, one now has to 
determine the actual value for pz, which depends on the 
parameters of the underlying noise model and should be 
chosen minimal (see appendix A), since this gives the 
strongest upper bound. The exact values for the bounds 
obtained in this way are however worse than the upper 
bound derived in Sec. IV ("I but as we will see in Sec. IV Fl 
the way of deriving the upper bound here will turn out 
to be well suited for all those cases where the initial state 
is only slightly entangled. 

Moreover, as it was the case for the lower bound, the 
derived upper bound on the lifetime of distillable entan- 
glement does only depend on the maximum degree of 
the graph and not necessarily on the number of particles 
N. We remark that the upper and lower bound on the 
lifetime of graphs states show a different dependence on 
the degree m of the graph. While the lower bound on 
the lifetime decreases with m, the upper bound on the 
lifetime increases with rn. We emphasize that this obser- 
vation applies only to the lower and upper bounds, and 
no definitive statement about the actual dependence of 
the lifetime of distillable entanglement on the degree of 
the graph can be made (although one may expect that 
the lifetime of entanglement decreases with the degree 
of the graph). The different dependences of the lower 
and upper bound can in part be understood by looking 
at the corresponding derivations. In particular, in the 
derivation of the upper bound t> the influence of both 
ax and ay noise is completely ignored. The influence of 
this kind of noise, however, strongly depends on the de- 
gree of the graph and is in fact responsible that e.g. the 
fragility of GHZ states depends on the number of parti- 
cles n4. That is, ax noise on all neighboring particles 
acts as az noise on a given vertex, and the noise accu- 
mulates. However, it is not straightforward to take also 
si'^'^ and Sy''^ in above analysis into account, as they lead 
to correlated noise when expressed in terms of az opera- 
tors (see Eq. (I44f) '). This implies that one could no longer 
consider separability properties of two-qubit maps inde- 
pendently but has to take correlations into account and 
thus consider a larger (or eventually the whole) system, 
thereby losing the main advantage of this approach on 
determining separability of the resulting state. 



E. Upper bound III: Partial transposition criterion 
for graph diagonal states 

In our third approach, we determine an upper bound 
on the lifetime by considering the partial transposition 
with respect to several partitions. Although this upper 
bound will be worse than the upper bound in Sec. IV CI 
(except for some singular cases), the ability to explicitly 
compute the partial transpose with respect to different 
partitions will enable us to compare the afore mentioned 
bounds with the exact critical values for the PPT crite- 
rion, at least for graphs with only few vertices {N < 10). 
In this case the techniques developed in this Sec. thus 
lead to stronger results for the lifetime of A^-party en- 
tanglement. For the following upper bound we make use 
of the fact that a A^-particle state is certainly no longer 
TV-party distillable if at least one of the partial transpo- 
sitions with respect to all possible bipartite partitions is 
positive. To this aim, we determine the eigenvalues of 
partial transposition of p{t) with respect to various par- 
titions. Since this is in general a rather complicated task, 
we will assume that decoherence of the particles is based 
on the same individual Pauli channel: 



3 

E 

J=0 



Pj{t)a'^pa!; 



(78) 



As it was already used in the derivation of the lower 
bound, under such a Pauli channel the graph state \G) 
evolves in time into a mixed state p = Yikev '^^''^ 1^) (^1 ' 
that is diagonal in the graph state basis \U)g'- 



P - 



E ^u\U)g{U\. 



(79) 



u<zv 



In the following we will make use of the following facts, 
whose proofs are postponed to the appendix B: 

(i) The diagonal elements Xu in Eq. l(7^ can be com- 
puted to be of the form 



Ac/ = Po ' 2^ 9l 



JV| ^ ^\U'\(TU' + U)\ JU'n(TU' + U)\ 

'''''xgf^'+^)\^'l, (80) 



where qi '■— for i = 1, 2, 3. In the case of the depolar- 
izing channel [q := qi = q2 = q3 = this simplifies 



to 



At,=pl^l ^ q\u'u(ru'+u)\ 



(81) 



U'CV 



We note, that in both expressions we have made use of 
the notational simplifications described in Sec. IV Al 

(ii) For any state p of the form H79() . i.e. that is 
diagonal in the basis \U)g according to some graph G, the 
partial transposition p'^^ with respect to some partition 
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A is again diagonal in the (same) graph state basis \U)g- 
In order to compute the corresponding eigenvalues, let 
r' = TaA" denote the adjacency matrix of the graph 
between the partition A and its complement A'^, i.e. 



^AA^ ^A^ 



r. 



Then: 



/^^■^ = E ^u\U}g{U\ with 



(82) 



(83) 



ucv 
Ikerr'l 



(_l)<^^^-> A 



(U+X+Y) , 



(x,y)e 
(kor r')^x(imr') 



where Ay S ^ is arbitrary with T'Ay = Y and the 
kernel ker or the orthocomplcment _L are taken with 
respect to the subspace ^^{A) spanned by the sets in A. 

(iii) In the case of small noise < = ^ < 1 for 
i = 1,2,3 the following estimation: 



qXu <: \u+k < - A[/ 

q 



(84) 



can be derived, where q = min((7i, (72, 93)- The same 
holds for Ajj+Affc and Ay+jv^+fc instead of Xu+k- 

Before coming to the upper bound let us give two ex- 
amples for formula H83|l: If T' is invertible, then kerF' = 
{0} and (kerr')-L = V{A) holds. Moreover ^ can be 
simplified by parameterizing ImF' with Y — V A2, where 
A2 C A: 

>^u-:^ E (-l)<^^^^^^V+^,+r'A.). (85) 



If j4 = {k} for a non-isolated vertex k £ V the eigenvalues 
of the partial transposition with respect to A are 

= 2 ^^'^ ^U+Nk + Xu+k — Xu+Nk+k) ■ (86) 

Similarly for the partial transposition with respect to the 
split A = {k,l} versus rest, where k,l G V are two non- 
adjacent vertices with linearly independent neighbor sets 
Nk and Ni, one obtains: 



\xeM+ 



u+x — E ^u+x 

X£M^ 



(87) 



where 



M+ = {0,k,l,k + l,Nk,Ni,Nk + Ni,k + Ni, 
I + Nk,k + I + Nk + Ni} and 

M- = {k + Nk,l + Ni,k + Nk + Ni,l + Nk + Ni, 
k + l + Nk,k + l + Ni} . 



If k and / are adjacent the same formula holds but with 
neighbor sets N'f^ = Nj. \ I and Nj^ = Ni \ k restricted to 
A'. 

Finally we note, that for GHZ diagonal states of the 
form H23() the positivity [non-positivity] of the partial 
transpose with respect to all possible partitions was al- 
ready a necessary and sufficient condition for iV-party 
separability [distillability] . For general graph diagonal 
states the corresponding PPT [NPT] criterion is only 
known to be a necessary condition for A^-party sepa- 
rability [distillability], whereas the sufficiency of these 
conditions is presently unknown. But for all partitions 
{A, A'^), for which the pure graph state |G) has Schmidt 
measure 1, i.e. it can be decomposed into the form 
\G) = ai|ai)^|6i)^' + a2|a2)^|&2)'^', the NPT criterion 
is also sufficient condition at least for the distillability 
of a {A, ^'=)-entangled state: If p'^^ in Eq. has a 
negative eigenvalue AJy, then the corresponding eigen- 
state \U)g has a Schmidt decomposition of the form 
\U)g = |G) = ai|a'i)^|5'i)^' + a2|a'2)^|6'2)^' and also 
a negative overlap g{U\p^'^\U)g < 0, which is sufficient 
for (^,A'=)-distillability 

These results can now be used to derive upper bounds 
to the distillable entanglement in graph states in the pres- 
ence of local noise described by a Pauli channel H18|l with 
Pi > for i = 1,2,3. For example, if one considers the 
split one- versus- rest, the eigenvalues X'jj of the partial 
transposition with respect to the corresponding partition 
A = {k} in Eq. H86(l can be bounded from below by 

X'u>{l + 2q-^)Xu. (88) 

Therefore the state p in H79|l is certainly PPT with re- 
spect to the partition A = {k} ii 1 + 2q ~ ^ > 0, i.e. 

i < g < 1. In the case of the depolarizing channel 
{q = 'i^iyi) this means, that no A — A'^— entangled state 
can be distilled from any graph state p if p{t) falls below 

^> = i 

For the partial transposition with respect to the par- 
tition A = {k,l} (see (|HZ|l), Eq. (|HH can be ap- 
plied twice (e.g. Xu+k+Ni > qXu+k > q^Xw ) in or- 
der to obtain estimations for the 'higher order' terms 
k + 1, k + Ni,k + l + Nk + Ni... in M+ and In this 

way one arrives at the condition 

A',7> (l + 4g + 5q2---4)V>0 (89) 
q q^ 

for the distillability of ^ — ^'^—entanglement. This 
means, that in the case of the depolarizing channel for 
q > 0.8457 or p < 0.0436 any graph state \G) will be- 
come PPT with respect to ^ = {k,l}- A closer com- 
parison with Sec. IV CI shows, that the upper bounds to 
A^-party distillable entanglement derived in this way are 
worse than the upper bound in Eq. (|60|l . For the one- 
versus-rest split, this can be understood for a general 
Pauli channel with > by rewriting the condition as 

min Pi > TT • (90) 

i=l,2,3 2 
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Then, due to q > 1, po must be the maximum in Eq. H60() 
and hence can only be larger than i if also ^ > j holds, 
which cannot be exceeded by the minimum in . Nev- 
ertheless, we think that the derivation can be of interest 
for other applications involving the partial transposition 
of graph diagonal states. In particular, it is an open 
question, whether for certain Pauli channels with pi > 
the conditions for PPT with respect to larger partitions 
might yield a stronger upper bound than the condition 
in Eq. (|6U|I . This will certainly depend on the solutions 
to the corresponding polynomial equation in q. But any 
upper bound derived with the use of (iii) will not depend 
on the topology of the underlying graph in question. By 
using a slight modification of the argumentation leading 
to the estimation in (iii) we will therefore discuss the 
example of the dephasing channel, for which a stronger 
upper bound can be provided, that conversely depends 
on the topology of the graph. In any case, the procedure 
to compute the eigenvalues of the partial transposition 
described in (iii) does not require the diagonalization of 
a 2^ X 2^-matrix and therefore allows the evaluation of 
the PPT criteria with respect to different partitions, as 
long as the vector consisting of the initial eigenvalues Xu 
(which is already of length 2^) is small enough to be 
stored and -in the case, that it occurs as a result of Pauli 
channel- as long as this vector can also be initialized fast 
enough. In order to illustrate the afore mentioned re- 
sults we have, for example, considered rings up to size 
= 10 suffering from decoherence due to the depolar- 
izing channel and examined the partial transpose with 
respect to all possible partitions. Fig. [T] depicts the crit- 
ical value for p, after which the state p first becomes 
PPT with respect to some partition, which implies that 
at this point the state p is certainly no longer iV-party 
distillable. For Fig.[7|the critical value Pcrit has also been 
computed, after which the state p has become PPT with 
respect to all partitions, i.e. after which p contains at 
most bound entanglement with respect to any partition. 
In contrast to the case of iV-party GHZ states, for which 
the one-versus-rest partition is the first to become PPT, 
the numerical results for small TV indicate that in rings 
this split seems to be most stable against decoherence 
due to noise described individual depolarizing channels 
and that the smallest eigenvalue of the partial transpo- 
sition with respect to these one-versus-rest splits {k} is 
given by XNk+k- 

In Fig. |S| we show representatives of the equivalence 
classes for connected graphs over N = 5, 6, 7 vertices dis- 
cussed in Q, that are most stable or instable, when ex- 
posed to noise described by individual depolarizing chan- 
nels. In this context we consider two graphs to belong to 
the same equivalence class if they can be transformed into 
each other by local unitaries and graph isomorphics. The 
latter corresponds to an exchange of particles, that maps 
neighboring particles onto neighboring particles. We note 
that in this special case of noise due to the same individ- 
ual depolarizing channel the notion of equivalence classes 
of graph states under local unitary transformations and 
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FIG. 7: For the case of particles in rings of size A'^ < 10, 
which individually decohere according to the same depolar- 
izing channel 11141 with parameter p: the critical value Pcrit, 
after which the first [last] partition becomes PPT A [□], the 
lower bound according to Sec. IV Bl and the upper bounds ac- 
cording to Sec. ED and Sec. EH 



N=5 



N=6 



N=7 




FIG. 8: Representatives of the equivalence classes under lo- 
cal unitaries and graph isomorphics of the connected graphs 
with A'' — 5,6,7 vertices f^. The first [second] column de- 
picts a representative of the class, that is the last [first] class 
of a given size A'^ to become PPT with respect to some parti- 
tion. Hence the first [second] column contains those graphs, 
for which the PPT-criterion indicates, that the A^-party dis- 
tillable entanglement contained in these states might be most 
stable [instable] (among all graphs with the same number of 
vertices). Similarly the third [fourth] column shows a graph 
of the equivalence class, that is the last [first] of a given size 
N to become PPT with respect to all partitions. Hence the 
third [fourth] column contains those graphs, for which the 
PPT-criterion indicates, that these states might be the last 
[first] to become A^-party separable. 



graph isomorphics (i.e. particle exchange) is meaning- 
ful, since the decoherence process itself is invariant under 
these operations. As it can be seen in Fig.|Slfor connected 
graphs on A^ = 5, 6, 7 the A^-party GHZ states seem to 
be the first that looses A^-party distillability. 

Finally we will consider the case of individual dephas- 
ing channels Vp = pp + [p + (JzPO'z) (i-e. po = 
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for which the estimation 



is no longer vahd in general. It is straightforward to 
see that 



holds for q — ^—^ < 1, since A 
^ i+P — ' 



pNq\U\_ 



(91) 

Similarly as 

in Eq. H88I) . we therefore can bound the eigenvalues of 
the partial transpose A'^; with respect to the partition k 
from above by 



A'c/>(l + gl^^-l+g-9-(l^^l+i))Ac 



(92) 



As depicted in Fig.|51 the above case of a ring {\Nk\ ==2) 
this inequality yields to the sufficient condition q > 
0.7549 [p < 0.1397] for all one-versus-rest splits to have 
PPT and hence yields a stronger criteria for iV-party dis- 
tillable entanglement than Eq. ()60|l . Note, that the lower 
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particles k and I have interacted, 
state \G) is thus given by 



\G) 



n u>^^\ 

{kd}eE 



The weighted graph 



(93) 



where the operations Uki are in this case given by 

[/,,^e-^(^<^'-'^^^0«(i<''-i"). (94) 

In contrast to this straightforward extension of the in- 
teraction picture for weighted graph states, no such gen- 
eralization of the stabilizer formalism (see Eq. (|38|l ) in 
terms of generators within the Pauli group is possible. 
In particular this implies that the results of Sec IV El are 
no longer applicable to weighted graph states. But in 
the following we will show, that the other techniques es- 
tablished in the previous sections to obtain lower and 
upper bounds on the lifetime of entanglement can in fact 
be extended to cover also weighted graph states. Ac- 
tually the following analysis will also hold for all states 
produced from Hki acting on an arbitrary product state 
\ip) = {ipi) . . . \4'n), which are not necessarily of the form 
lip) = \+)'^^ . Nevertheless, for the sake of simplicity we 
will restrict the following to this case. 



Lower bound on lifetime 



FIG. 9: For the case of particles in rings of size A'^ < 10, 
which individually decohere according to the same dephasing 
channel with parameter p: the critical value Pcrit, after which 
the first [last] partition becomes PPT A [□] , the lower bound 
according to Sec. I V Bl and the upper bounds according to 
Sec. ED and Sec. EH 

bound actually coincides with the computed critical val- 
ues for p, after which the ring first becomes PPT with 
respect to some partition. 

F. Generalization to weighted graph states 

In this section we extend the previous results to a more 
general class of initial states, the so called weighted graph 
states. The graph states discussed so far arise from the 
Ising type interaction Hamiltonian i/fej = |1)'^(1| (g) |1)'(1| 
(see Ea. (|36f) ') acting on a collection of particles V in the 
ax eigenstate for a fixed time (p = ir according 

to some interaction pattern specified by the graph. We 
will now allow the particles to interact according to the 
same Hamiltonian Hki but for different interaction times 
ifki- This corresponds to the situation of a disordered 
system as it occurs e.g. in a spin glass or semi-quantal 
Boltzmann gas. The interaction pattern can similarly be 
summarized by a weighted graph, in which every edge is 
specified by a phase ipki corresponding to the time the 



In order to obtain a lower bound on the lifetime of 
weighted graph states, we again provide an explicit pro- 
tocol which allows the distillation of maximally entangled 
states between all neighboring pairs of particles and thus 
to create arbitrary Ai"-particle entangled states. In fact, 
we make use of the same protocol as in Sec lVB| however 
the analysis turns out to be different. To be specific, we 
consider the state p(t) which is obtained from a weighted 
graph state \G) subjected to decoherence -described by 
individual Pauli channels- for time t. We perform mea- 
surement in the eigenbasis of on all but particles 
k, I and determine the condition when the resulting re- 
duced density operator pki is distillable. We denote by 
Pq = |0)^(0|,Pf = projectors acting on parti- 

cle j and by Po,Pi the corresponding superoperators, 
i.e. Pqp — PqpPq. Similarly, we denote UkiQ = UkigUli 
the superoperator corresponding to the unitary opera- 
tion Uki- Note that the entanglement properties of the 
resulting state do not depend on the specific measure- 
ment outcomes. For notational convenience we restrict 
our analysis to the case, where the measurement result 
is obtained on all measured particles. Taking noise de- 
scribed by some Pauli channel Dp — ^^^QPiaipcTi, we 
thus have to consider the (unnormalized) state 

VkVi n P^T^i n Uab\+)''{+l (95) 

j^k,l {a,b}eE 
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Using that Pqo-q.s = o-q ^Pq and PoO'i,2 = o'i,2Pi we can 
rewrite P^'DjQ and obtain 

P^V, Q = (MiPa + M{Pi) Q, (96) 

where 

MiQ = piaiQdi + p2(J2Qi^2- (97) 

Choosing the computational basis |[/)^\{'=''} = 
(y^\^)^^^^'^^ on the measured particles, we thus can write 

\{P^V,g^ E W^uPlj^S (98) 

j^kj UCV\{kd}j^kd 

The projector commutes with the unitary operations Uat 
and we therefore obtain that Eq. H95|l can be rewritten 
as 

v,vi u^'u. n ^'^^ 

UCV\{kd} 3^k,l {a,fc}eE 

|++)^-'(+ + |® |C/)^\^'=^'>(f/| . (99) 

Note that Uab leaves |[/)^\{'=''l invariant and it is thus 
sufficient to consider only Uab that act on particles fc, I 
and/or their neighbors, i.e. the set / = Nk U NiU kU I. 
For all other particles j ^ I we thus have expressions 
of the form MujlUj)^ {Uj\ — |0)-'(0|, i.e. these particles 
factor out. It follows that the reduced density operator 
Pfci which is obtained by tracing out all but particles fc, I 
only depends on particles in the set / but not on the 
other particles j € / or errors (noise) effecting these other 
particles. This already shows that the lower bound on 
distillability for weighed graph states only depends on the 
(degree of the) corresponding interaction graph as well as 
the weights of the edges, but is independent of the size 
of the system N as long as the degree of the graph itself 
does not depend on TV. In particular, only the subgraph 
of particles j € I determines the entanglement properties 
of the reduced density operator pki- 
We have that pki is given by 

ucv\{kj}jei 

I ++)'■'(+ + I ®|C/)^\^^-''>([/|) , (100) 

where the partial trace has to be performed for the re- 
maining neighboring particles I oik and / only. Thus the 
effect of noise can be localized to the region / around the 
edge {fc, V\ in question. In principle, one can now obtain 
the explicit form of p^i for a given (weighted) graph and 
determine the condition for pi until when the reduced 
density operator pki has non-positive partial transposi- 
tion and remains thus distillable. The explicit formula is 
however rather complicated and not particularly illumi- 
nating. For the example of a depolarizing channel, it is 



clear that for smaller values of Lpu (i.e. a weaker edge 
between particles k and 1) one obtains stronger threshold 
values on the parameter p than given by Eq. (g^, i.e. a 
shorter lifetime. 

What is however more important in our context is that 
also for weighted graph states the lower bound on lifetime 
of distillable entanglement only depends on the (degree) 
of the corresponding interaction graph, but not on the 
size of the system N . Although the actual values of the 
lifetime will depend on the specific weights of the edges, 
for cluster-like and similar graph states there will be no 
scaling behavior with N . Moreover, in many cases such 
as rings, the edge with the smallest weight will give rise 
to the strongest threshold value condition and will thus 
determine the lower bound on the lifetime of distillable 
iV-party entanglement. Actually, it is sufficient if one 
can create maximally entangled pairs between pairs of 
particles in such a way that there exists a path between 
each pair of particles (i.e. entanglement between all pairs 
{fc, 1} where the edges {fc, 1} form a maximally connected 
graph). Thus the state is already A^-party distillable, if 
the subgraph is connected, that is generated by all those 
edges, from which a Bell pair can be distilled. This im- 
plies that some edges in the original graph -even if they 
are very weak- may not play a role if there exists an- 
other way to obtain singlets between all relevant pairs. 
For instance, if one considers a graph corresponding to 
a lZ?-cluster state, where each edge has weight tt, and 
one adds in addition an edge 1, k with small weight (pi^, 
then it is not necessary that entanglement between par- 
ticles 1, k can be distilled (although the two particles are 
neighboring ones according to the graph G), but it is 
sufhcicnt to distill entanglement between all pairs of par- 
ticles fc. fc + 1. 



2. Upper bound on lifetime 

The first method to obtain an upper bound to the life- 
time of entanglement certainly also holds for arbitrary 
weighted graph states, since it is independent of the ini- 
tial state and reflects the time after which the decoher- 
ence process itself has become entanglement breaking. 
Conversely the upper bounds derived in this way can- 
not take into account whether the initial state is only 
slightly entangled or not. In Sec. IV Dl it turned out that 
for ordinary graph states the upper bound is weaker than 
the first one derived in Sec. IV CI In contrast we will 
show in the following that the upper bound presented in 
Sec lVDI will give tighter upper bounds to the entangle- 
ment in weighted graph states, which contain vertices k 
with only small interaction phases ^pki to all their neigh- 
bors I G Nk- We use again the dynamical description 
of the weighted graph state \G) given by Eq. (|^ . The 
influence of (phase) noise on the entanglement properties 
of Uki can be determined in a similar way as in Sec lVDl 
where we had a fixed angle ipki = tt for all {k, 1} G E. 
We remark that the upper bound obtained in Sec IV Dl 
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for a general graph state is also valid for all graphs of 
the same kind where the edges are weighted. This is due 
to the fact that the operations Uki are most resistant to 
noise (i.e. remains entangling), if the angle is fki = i", 
because in this case the operation is -in the ideal case- 
capable to create maximally entangled states, while for 
all other values of ipki only partially entangled states can 
be created. 

This observation immediately leads a way to obtain 
stronger upper bounds on the lifetime for weighted 
graph states. To this aim, one determines the thresh- 
old values when the state p{t) (Eq. be- 
comes separable. The value of p^' now depends on 
(fiki- One finds that the corresponding density opera- 
tor Dk^kMAiPzY^^^'KipfY^^^'^) in Eq. m again has 
support only in a four dimensional subspace and is given 

by 



D 



1 

kik2iii2 = X! 

hj=0 



(101) 



with 



Aoi - \ {l + {pf)^) (l-(pf)^ 
l{l-{pf)T^){l + {pf)T^ 
\{^-{pf)^){i~{pf)^ 



Aio 
All 



(102) 



The orthogonal states are given by {ai'^^^'Y 

(ct^^)^>oo), where 



|$oo> = |6)fe| 



|0)fc|i), + |i)fc|0),+e^'^-|i)fc|i)(103) 



with |6)fc = |00)feifc2, = 1 11) and similar for particles 
lil2- We have again, that Dk-^k2iii2 separable if and 
only if the partial transposition is positive, which leads to 
the threshold value . The separability of the weighted 
graph state can then be determined in a similar way as in 
Sec lVDI In order to make the operation Um separable, 
az noise with (p^')i/l^fcl and (p^')^/l^'l is required at 
vertices k and I. At a given vertex ko, this leads to a 
required total value of p^" = mini^N^^ p^f' such that all 
operations Ukgi become separable. The threshold value 
p>, below which the state is fully separable, is finally 
obtained by taking the minimum over all (that is over 
all vertices). That is 



P> 



mmp^, 

k 



(104) 



and the state p is certainly separable for pz < p>- For 
the different decoherence models £ allowing for an ex- 
traction of a dephasing component (see Eq. H125() ') one 
finally has to insert the relation between the dephasing 
parameter pz and the noise parameter of £ in order to 
obtain the announced upper bound on the lifetime for 
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FIG. 10: For the case of individual depolarizing channels with 
parameter p the threshold value Pcru = 2+p~ '^l' lOSJ) 
as a function of the interaction phase ipki £ [0, n] for edges 
{k,l} between two vertices k and I with the same increasing 
degree m = \Nt\ = lA^il = 2, 3, . . . , 10. The horizontal line 
depicts the upper bound according to Sec. I V Ul in this case. 



weighted graph states. Fig. [TUI depicts the critical value 
for the depolarizing parameter p in Eq. I|14|l as a func- 
tion of the weight ifki at the edge in question. As it 
was already mentioned in Sec. IV Dl for a fixed phase (pki 
the obtained upper bound on the lifetime decreases with 
the degree of the neighboring particles k and I (in con- 
trast to the corresponding behavior of the lower bound 
in Sec. IV B| ). Moreover Fig.llOlshows. that for any degree 
m there always exists a range of values for ipki , for which 
the analysis of this Sec. provides a stronger upper bound 
Pcrit than the condition derived in Sec. IV CI 



VI. BLOCK WISE ENTANGLEMENT AND 
RE-SCALING 



Entanglement is a concept which can only be defined 
between subsystems of the whole system. In our previ- 
ous analysis, we have identified subsystems with parties, 
i.e. we have investigated the lifetime of true A'^-party 
entanglement. One can however also consider a slightly 
more general concept, where subsystems are formed by 
a collection of several parties (see Sec. IIII|l . Also in this 
case one can investigate entanglement properties between 
M such subsystems. That is, one can consider a parti- 
tioning of the A'^-party system into M < N groups and 
investigate the (distillable) entanglement between these 
M groups, where each of these groups consists of one or 
several of the initial parties. Whenever M < N, one can 
have that the state is still M-party entangled although 
it contains no longer A'^ party entanglement. Consider- 
ing such coarser partitions allows one to investigate the 
change of the kind of entanglement in time and to deter- 
mine an "effective size" of the entanglement present in 
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the system. One can determine for each kind of entan- 
glement the corresponding hfetime. 

A. Block wise entanglement: Distillability and 
lifetime 

Given an N party system, we consider a partitioning 
of the N parties into M < N groups (M-partitioning) . 
Parties within a given group are allowed to perform joint 
operations are considered as a single subsystem with a 
higher dimensional state space. We are interested in the 
entanglement between these M subsystems. A density 
operator p is called M-party distillable with respect to 
a certain M-partitioning if from (asymptotically) many 
copies of p one can create some irreducible entangled 
pure state by means of local (in the sense of the M~ 
partitioning in question) operations and classical commu- 
nication. Similarly, a density operator p is separable with 
respect to a certain M-partitioning if it can be written 
as convex combination of products states (in the sense of 
the M-partitioning in question). 




FIG. 11: For blockwise entanglement different partitionings 
of particles into groups are considered. 

We say that a density operator is M-party distill- 
able if it is distillable with respect to at least one M- 
partitioning. It is obvious that if p is M-party distill- 
able, it as also M'-party distillable for M' < M. The 
maximal possible M such that a density operator p is 
M-party distillable can be interpreted as a measure of 
the "size" of entanglement, as it provides the maximal 
number of subsystems (groups, blocks) which are (dis- 
tillable) entangled. We will investigate the lifetime of 
M-party distillable entanglement for all M. We remark 
that M-partitionings might be completely different in 
nature. Consider for instance a system consisting of 90 
parties Ak,k = 1,...,90. Possible S-partitionings in- 
clude for instance [Ai] — {A2) — (A3, A4, . . . , Ago) as well 
as (Ai, . . .,^30) - (A31, . . . , Ago) - (Aei, . . . , Ago). For a 
detailed description of M-partitionings as well as neces- 
sary and sufficient conditions for M-party distillability 
we refer the reader to Ref. ,11.] . 

Equivalently one may say that we are considering a 
coarse- graining of the partition and investigate the entan- 
glement properties under coarse graining. As a particular 
instance of such coarse graining we will investigate par- 
titionings that correspond to re-scaling of the size of the 



subsystem as it is used in statistical mechanics. Consider 
for instance N particles which are arranged on a regular 
rectangular (two-dimensional) lattice. The finest parti- 
tion corresponds to considering each particle individually 
as a single subsystem. Coarsening of the partition may 
e.g. take place by considering blocks of n x n particles 
(arranged as a square) as a single subsystems. That is, 
for n = 2 one considers a specific M = N/A partition- 
ing, while for arbitrary n we have M = N/'n? groups of 
parties/subsystems, each formed by v? particles. This 
concept is also illustrated in Fig. We will investigate 
how (distillable) entanglement changes under such re- 
scaling and determine asymptotic properties of the life- 
time of M-party entanglement for macroscopic number 
of particles, — > 00. 



B. Lifetime of GHZ states under coarse graining 

1. Lifetime of entanglement in large T limit of reservoir 

As in Sec. IIVI we consider the lifetime of distillable 
entanglement for GHZ states when each particle is in- 
dividually coupled to a thermal reservoir with T — > 00, 
described by a depolarizing quantum channel. In order to 
determine the lifetime of M-party entanglement, we can 
make use of the results obtained in Sec. IIVI together with 
the classification of GHZ-diagonal states of Ref. (llj . 
Recall that the partial transposition with respect to a 
group (subsystem) Bk which contains exactly k parties 
is positive, p[t)^^k > if and only if < 2\k (see Eqs. 
(|25l26|l ). In addition, we have Ai > A2 > . . . > X[n/2] 
(see Eq. l|?7j) '). This implies that the size of the subsys- 
tem, i.e. the number of particles that are contained in 
a subsystem, determine when the corresponding partial 
transposition becomes positive. We have that the par- 
tial transposition with respect to a single party, p{t)^^>' , 
is the first one that becomes positive, while the partial 
transposition with respect to a larger subsystem is more 
stable, i.e. becomes positive at a later time. 

For a given M-partitioning, above observation im- 
mediately allows one to identify the bipartite partition 
(which contains the M-partition in question) which de- 
termines the lifetime of distillable Af -party entanglement. 
In particular, the partial transposition with respect the 
subsystem that contains the smallest number of parties 
is the first one to become positive. We use that a nec- 
essary -and in the case of states of the form p{t) we 
deal with in this case also sufhcient- condition for M- 
party distillability is that the partial transposition with 
respect to all subsystems Bk forming the M-partitioning 
is non-positive. To be exact, one also needs that all 
partial transposition with respect to groups formed by 
several such subsystems Bk are also non-positive, which 
is in our case however automatically fulfilled as partial 
transpositions with respect to larger groups are more sta- 
ble than with respect to smaller groups. If only one of 
all these partial transpositions with respect to various 
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subsystems is positive, it follows that the state p{t) is 
no longer Af-party distillable (see Ref. 0). Thus we 
have that the lifetime of M-party distillable entangle- 
ment with respect to a given M-partitioning is deter- 
mined by the size of the smallest subsystem of the corre- 
sponding M-partitioning. For instance, if one considers 
an arbitrary M-partitioning {M < N) that contains as 
one subsystem a single party, say Ai, we have that the 
lifetime of entanglement with respect to this partition- 
ing is completely determined by the partial transposition 
with respect to party Ai, i.e the condition < 2Ai. In 
particular, the lifetime of distillable Af -party entangle- 
ment with respect to any such partitioning is exactly the 
same as for iV-party entanglement (corresponding to a 
iV-partitioning where each subsystem contains a single 
party). 

It follows that for a given M, M-party entanglement 
with respect to a specific M-partitioning has longest life- 
time if all groups have (approximately) the same size. For 
a minimal group size of m particles, a A'^-particle GHZ 
state can contain at most M = [N/m] such groups of 
size m. This allows one to obtain the maximum life- 
time of M-party entanglement which is determined by 
< 2Xm (Eq. ^ with k = m). 

Determining the threshold value pcrit involves the so- 
lution of a polynomial equation of degree N, which can 
be done numerically in an efficient way. One can how- 
ever also determine analytic lower and upper bounds on 
the lifetime of M-party entanglement. One obtains an 
upper bound on the lifetime of M-party entanglement if 
one approximates Am by some Xm < and investigate 
the condition 



(105) 



as in this case automatically also p^ < 2Xk and thus the 
partial transposition with respect to a group that con- 
tains m parties is certainly positive. We can e.g. choose 



Am = (l-p)™(l+p)^-™/2^+\ 



(106) 



which obviously fulfills Xk < Xk- The condition Eq. IjlOSf) 
can then be rewritten as 



m< N 



\og[2p/{l+p)] 
logiil -p)/ {I +p)Y 



(107) 



Recall that in the (optimal) case where all subsystems 
have the same size we have M = [N/m]. Thus we find 
that a A'^-party system is certainly no longer M-party 
entangled if 



M > 



log(l ~p) ~ log(l -t-p) 
log(2p)-log(l+p) 



(108) 



Eq. H108|l is a central result in our analysis of the prop- 
erties of the lifetime for GHZ states under re-scaling. We 
have illustrated Eq. ((TU5|l in Fig. □ and H 

On the one hand, Eq. (|108|l provides an upper bound 
on the lifetime ktm of M-party entanglement in the sys- 
tem. This upper bound can be obtained by determining 



ktm for a fixed M from Eq. H108|) , where Kt is again given 
hy p = e""*. On the other hand, for a fixed time t one 
can determine from Eq. (|l()8|l the maximum M of distil- 
lable multipartite entanglement in the system. That is, 
the effective size of entanglement after a time t can be ob- 
tained this way. One observes (see Fig.^andlSl) that the 
maximum M rapidly decreases with t. For small times, 
i.e. nt <C 1, one finds that M scales as 



M sa -21og(Kt)/(Ki). 



(109) 



For Kt > 0.8049 (which is obtained from Eq. H108|) by 
setting M — 2), we have that also 2-party entanglement 
disappears. In fact, the state becomes fully separable as 
all partial transposes are positive (which is a sufficient 
condition for separability for such states EJ). We em- 
phasize that the upper bound on the lifetime of M-party 
entanglement is independent of the number of particles 
N, in particular it is also valid for a macroscopic number 
of particles and even for N oo. 

On the one hand this implies that even if the size of 
the groups m = N/M goes to oo (for a total number 
of particles — s- cx)), the maximum number of groups 
M that can remain entangled after a time t is finite^ 
i.e. the maximum effective size of entanglement is lim- 
ited for any time t. On the other hand it follows that in 
the limit N oo any partitioning in groups with finite 
size m leads to a vanishing lifetime of the corresponding 
M — N/m party entanglement (see Eq. HlOTfl ). Only if 
one considers the limit where the number of subsystems 
M is fixed as iV ^ cxd, i.e. the size of each of the groups 
itself is macroscopic and m — > oo, one obtains that the 
lifetime of the corresponding M-party entanglement is 
finite. We remark that above results also enable one to 
obtain (upper bounds) on the lifetime of iV-party entan- 
glement by considering m = 1 and M — N (see also 
SecHg. 

In an analogous way one can derive a lower hound on 
the lifetime of M-party entanglement. To this aim, one 
uses that if the partial transposition with respect to the 
smallest subsystem is still non-positive, also all other rel- 
evant partial transpositions of the system with respect 
to all other subsystems (and combinations thereof) are 
non-positive. For states of the form Eq. I|23ll this ensures 
that the state p{t) is M-party distillable 11]. Thus we 
can derive an analytic lower bound on the lifetime of M- 
party distillability by considering the condition Eq. I|26|) 
with k — m and upper bound A„i by some A^- We have 
that if p^ > 2 Am > 2 Am, then p{t) is certainly M-party 
distillable with respect to a partition which consists of 
M subsystems of size m = [N/M] . We choose 



A:„^2(1-p)"(1+p)^-™/2 



■)N+1 



(110) 



which can readily be checked to fulfill Am > Am- The 
condition p^ > 2Am can be rewritten and one finds that 
for 



M < 



log[2(l-p)]-log(l+p) 
log(2p)-log(l+p) 



(111) 
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all partial transposition with respect to this M- 
partitioning are certainly non-positive, which already en- 
sures that the state p{t) is M-party distillable. 

At first sight, our results seem to contradict the ones 
obtained by Simon and Kempe (181] . They observed that 
the threshold value for p when considering the partition 
{N/2) — {N/2) decreases with the size of the system N. 
Based on this observation, they conclude that GHZ states 
of more particles are more stable against local decoher- 
ence. However, as pointed out in the discussion above, 
the effective number of subsystem that remain entangled 
decreases with time. The entanglement becomes bipar- 
tite when approaching the threshold value found by Si- 
mon and Kempe. In fact, the lifetime of genuine (distil- 
lable) A^-party entanglement decreases with the size of 
the system N. 



2. Quantum optical channel 

One can now perform a similar analysis of the lifetime 
of M-party entanglement for more general couplings of 
the individual particles to the environment described by 
a general quantum optical master equation. We have 
already determined the condition such that the partial 
transposition with respect to a group of k parties is pos- 
itive in Scc lIVBI (see Eq. (j^ ). Following the hnc of 
argumentation if Sec I VI B"T1 and using the notation of 
Sec lIVBl it is straightforward to obtain an upper bound 
on the lifetime of Af -party entanglement. In particu- 
lar, one has that positivity of at least one of the partial 
transpositions with respect to a specific subsystem en- 
sures that p{t) is no longer M-party distillable. That 
is, when considering (as in the previous section) a M- 
partitioning of the A'^-party system into M subsystems, 
each of size m = N /M , an upper bound on the lifetime is 
given by Eq. ^A'A\ with k = m. In fact, also in this case we 
have that the subsystems that contains the smallest num- 
ber of parties gives rise to the strongest condition on the 
lifetime of M-party entanglement, i.e. the corresponding 
partial transposition is the first one to become positive. 
This can be seen by considering the condition for posi- 
tivity of the partial transposition with respect to a group 
of k particles given by Eq. [h^ /2f <\k^N^k- We 

have that for 1 < fc < [N/2] 



from Eq. (|33|l that p{t) is certainly no longer M-party 
distillable entangled if 



^k^N~k > ^k+l^N-k- 



(112) 



(which can be checked by direct computation and consid- 
ering separately the cases (ac) > and (ac) < 0). From 
this observation the claim already follows, as Eq. H112() 
ensures that if p{t) has positive partial transposition with 
respect to 1 parties, it automatically has also positive 
partial transposition with respect to k parties. Hence 
the subsystem with the smallest number of parties de- 
termines the lifetime of M-party entanglement. Using 
again that for a fixed M one obtains the longest lifetime 
if all M groups have the same size m = [N/M] , it follows 



M > 



log(ac)-log[(l-a)(l-c)] 
log(ac) + Bt 



(113) 



where a, 5, c are defined by Eq. H3U|) . Thus an analogous 
discussion as in Sec lVIB Tl apphes, i.e. the scaling 
properties of the lifetime of GHZ states with the size 
of the system A^ and the number of subsystems M is 
similar for different couplings of the particles to the 
environment. Only in singular cases (such as zero 
temperature corresponding tos = Oors = 1), one 
discovers a different behavior. 

Considering general graph states instead of GHZ states 
a detailed analysis of their blockwise entanglement has 
not been accomplished yet. Nevertheless, the scaling be- 
havior of Af -party entanglement is restricted to a range 
between the upper and lower bounds of Sec. EI which in 
the case of cluster and similar graph states were shown 
to be independent of the number A^ of particles. In this 
sense also the scaling behavior of blockwise entanglement 
in these states must be essentially independent of the size 
of the system. 



C. Lifetime of encoded entangled states 

Until now we have seen two different kind of scaling 
behaviors with respect to the number of particles A^. 
For GHZ states, we found that the lifetime of distillable 
entanglement decreases with A^, while for cluster states 
(and similar graph states) we have that the lifetime is in- 
dependent of A^. In this Section we will show that certain 
states can show a different scaling behavior, namely that 
the lifetime of (block-wise) entanglement can even in- 
crease with the number of particles A^. Examples of such 
states are provided by encoded entangled states, i.e. en- 
tangled states which are formed by logical qubits where 
each of the logical qubits corresponds to the codewords 
of a (concatenated) quantum error correction code. We 
find that entanglement between M such logical qubits is 
maintained. Each of the logical qubits forms a subsys- 
tem of size m, and we consider entanglement between 
M such subsystems. We will show that as m increases 
(but is still of finite size), there exist states such that the 
maximum number of subsystems M that remain entan- 
gled increases. In addition, for a fixed M we have that 
the lifetime of M-party entanglement increases as the 
block-size m increases and tends to infinity as m — s- cx). 
This shows that (encoded) macroscopic entangled states 
-even of GHZ-type- can persist for long times. Note that 
this behavior is in contrast to non-encoded GHZ states, 
were on the one hand the maximum number of subsys- 
tems M that remain entangled after a certain time is 
finite and on the other hand also the lifetime of M party 
entanglement is finite (even as m — > oo). 
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One can interpret these results in the sense that the 
time in the encoded system is slowed down as compared 
to the time in the original system. This provides an alter- 
native view on quantum error correction and allows one 
to understand why encoded macroscopic superposition 
states can be produced and maintained on a quantum 
computer. 

1. Quantum error correcting codes 

We consider two orthogonal states of m qubits, 
|0l),|1l) € ((y^)®™ which correspond to codewords of 
some error correcting code and constitute the basis state 
of a "logical qubit" . In the following we will consider 
an optimal error correcting code which allows to correct 
an arbitrary error on one of the particles and uses five 
physical qubits to encode one logical qubit, e.g. the five 
qubit Steane code In the following discussion we 

will assume that each of the physical qubits is coupled 
to an independent environment and the the individual 
coupling is described by depolarizing quantum channels 
T>{p) (corresponding to the coupling to a heat bath in the 
large T limit), Eq. fl^ . with p = e"*^'. We consider the 
evolution of an arbitrary state of a single logical qubit 

\^L)^a\0L)+f3\lL), (114) 

under the influence of depolarizing channels acting on 
each of the individual particles, i.e. 

g{t)=Vi{p)V2{p)V3{p)V^{p)V,{pML){M- (115) 

The action of the map Vk on physical qubit k is such 
that with probability 

g=(3p+l)/4, (116) 

no error occurred, while with probability (1 — q) the 
qubit was effected by some error. In particular, we have 
that one of the three possible errors described by ak, 
k = 1, 2, 3 occurred with probability (1 — 5) /3 = (1— p)/4. 
Considering now the logical qubit consisting of five phys- 
ical qubits, we know that there exists a sequence of oper- 
ations (error syndrome measurement followed by a cor- 
rection step depending of measurement outcome) such 
that the state of the five qubits remains in the subspace 
spanned by {|0l), |1l)} and the logical qubit remains in 
the initial state \4>l) as long as no or only a single er- 
ror in one of the physical qubits occurred. That is, with 
probability 

qL = q' + 5q\l-q), (117) 

no or only a single (correctable) error happened, while 
with probability (1 — q^) the logical qubit was effected 
by some error. By applying (correlated) random unitary 
operations on the subspace spanned by {|Oi),|li;,)} at 
t = and t, one can achieve that the errors can again be 



described by white noise acting on the logical qubit, i.e. 
a map of the form 

1 - ^ 

T^L{pL)pL^PLPL + ^^^ai'^^pA''\ (118) 

fc=0 

where pl ~ {'^<iL ^ l)/3 and cr[.'^'' denote Pauli opera- 
tors acting on the logical qubits, e.g. tij^^ = |0l)(1l| + 
|1l)(0l|- The parameter is related to the initial p via 

PL = (3j5+ 1)4(4 -3p)/192- 1/3. (119) 

That is, the action of the decoherence process on the log- 
ical qubit can (after performing a correction step plus 
depolarization) be described by a depolarizing channel 
acting on the logical qubit, where the parameter p^ can 
be obtained from p -the parameter describing the de- 
coherence process of the individual physical qubits- by 
Eq. (|119|l . We thus have that a logical qubit where 
each of the particles is subjected to decoherence for time 
Kt = — log(p) behaves as if it was subjected to decoher- 
ence described by a depolarizing channel acting on the 
logical qubit for time kIl = — log(pL). We remark that 
if p > 0.82517 [nt < 0.1921658), we have that pl > p 
{kIl < nt), that is the logical qubit is less effected by 
decoherence as a physical one. In other words, the effec- 
tive time tL for which the decoherence acts on the logical 
qubit is smaller than the physical time and thus the de- 
coherence process on the logical qubit is slowed down. 

In a similar way, one can consider logical qubits which 
are formed by codewords of concatenated quantum error 
correction codes with k concatenation levels. A logical 
qubit consists in this case of 5*^ physical qubits. This 
follows from the fact that at concatenation level j, each 
logical qubit at level j — 1 is replaced by five such logical 
qubits (of level j — 1) which form the new logical qubit at 
concatenation level j . Following the same reasoning as in 
the case of a single concatenation level, one obtains that 
a logical qubit at concatenation level j is subjected to 
decoherence described by a depolarizing quantum chan- 
nel acting at the logical qubit, where the noise parameter 
qj is related to the noise parameter qj-i of the (logical) 
qubits at concatenation level j — 1 by 

9j =91-1+59^4-1(1 -g.,--!), (120) 

where qj = {3pj + l)/4 and qo — q {po — p) corresponds 
to decoherence acting on the physical qubits. That is, at 
each concatenation level the effective time tj for which 
the decoherence acts on the logical qubit of level j de- 
creases as compared to the physical time t = to as long 
as Kt < 0.1921658. For nto = 4/3eo < 1 we have that 

e, « 10e2_i, (121) 

from which follows 

Kt, « ^ . (122) 
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That is, the effective time ntj for the decoherence acting 
on the logical qubit is drastically decreased. For instance, 
if Kto ^ Kt = 0.01, we have that Kti w 7.5 x 10~^, Kt2 ~ 
4.22 X 10'^ Kt3 « 1.33 X 10-^°, nU ~ 1.34 x lO'^^, nt^ « 
1.34 X 10-3"^, Kig ~ 1.35 X 10-'^3. Recall that the num- 
ber of physical qubits that form a logical qubit at con- 
catenation level j is given by 5^, which corresponds 
to 5, 25, 125, 625, 3125, 15625 respectively. In particular, 
this implies that by using 625 particles to encode one 
qubit (i.e. j = 4), the effective time as compared to 
the physical time t can be decreased by a factor of 10^^. 

2. Lifetime of block-wise entanglement of encoded states 

Using the effective change in the timescale for logical 
qubits, it is now straightforward to determine lifetime of 
encoded entangled states. We consider block-wise en- 
tanglement, i.e. entanglement between M subsystems 
where each subsystem consists of m qubits. We consider 
lifetime of different kinds of entangled states which are 
formed by logical qubits and where each physical qubit 
is subjected to decoherence described by a depolarizing 
quantum channel (Eq. H14|l ) as in Sec lIV Al We empha- 
size that a logical qubit behaves in this context in exactly 
the same way as a physical qubit, the only difference is 
that the effective time (or noise parameter pl) is dif- 
ferent. This implies that the results obtained for GHZ 
and graph states can be directly used to obtain lifetime 
of such encoded entangled states. 

Consider an entangled state of GHZ-type which is 
given by 

\GHZl) = ^(|0l)«^^ + IIl)""^'), (123) 

where {|0l), |1l)} are codewords of (concatenated) quan- 
tum codes (with j concatenation levels) formed by m 
physical qubits which represent a logical qubit. When 
considering block-wise entanglement between M blocks 
of m = 5^ qubits, this state behaves in exactly the 
same way as a M particle GHZ states (Eq. H22|l ) con- 
sisting of physical qubits where one considers blocks of 
size I. The only difference is that the effective timescale 
tj (effective noise parameter pj ) is changed according to 

Eqs. dnoini. 

When considering the original GHZ state and Kt — 
10~^, we obtain from Eq. IjlOSI) (i) After a time nt — 
10"^, the maximum number of blocks that can be entan- 
gled is upper bounded by M = 1057; (i) For M = 1057, 
Kt = 10~^ provides an upper bound on the lifetime of M 
blocks (where in both cases each block may have arbi- 
trary size, i.e. m —^ oo). 

For an encoded GHZ state (Eq. H123|) we have the fol- 
lowing results: (i) The maximum number of blocks of size 
m = 5^ (i.e. logical qubits) that is entangled after a time 
Kt — 10^^ is determined by Eq. (|108|) . where p in this 
equation is given by e""*^ and Ktj is the effective time. 
We find that for j = 1,M = 2.103 x 10'^; j = 2, M = 



6.195 X lQ6;j = 3, M = 3.510 x 10", i.e. the number 
of blocks of fixed size m = 5^ that remain entangled af- 
ter some time Kt is drastically increased, (ii) An upper 
bound on the lifetime of M = 1057 blocks consisting of 
logical qubits of size m = 5^ is provided by Ktj = 10^^, 
where Ktj is the effective time. One can determine the 
time Kt (which provides an upper bound on the lifetime 
of such systems) using the recursive formula Eq. (|120|l . 
One finds that for j = 1, Kt = 0.0382; j = 2,Kt = 0.0778; 
j ^ 3,Kt ^ 0.1149; j ^ A,Kt ^ 0.1431; j = 5,Kt 
0.1621. Again, one sees that the lifetime of M-party en- 
tanglement is increased as compared to the original GHZ 
state, although here we considered only blocks of finite 
size m = 5-' , while we allowed for blocks of arbitrary size 
TO — > cx) in the case of (original) GHZ state. We remark 
that one can only expect that the encoded system has 
longer lifetime as compared to the original GHZ state as 
long as Kt < 0.1921658, since only in this case tj < t. 

In a similar way, the lifetime of M-party cluster- and 
graph states (formed by logical qubits) is enhanced when 
considering such states formed by logical qubits. 



VII. SUMMARY AND CONCLUSIONS 

In this paper we have investigated the lifetime of (dis- 
tillablc) entanglement under the influence of decoherence. 
We found that the qualitative behavior of different kinds 
of entangled states are largely independent of the specific 
decoherence model. In particular, we found for (essen- 
tially) all decoherence models with individual coupling 
of particles to (independent) environments that the life- 
time of GHZ states decreases with the size of the sys- 
tem. On the other hand the lifetime of cluster states 
and graph states with a constant degree (which does not 
depend on N) is independent of the number of particles 
N. The last observation can even be extended to all de- 
coherence models which correspond to some correlated 
but localized noise, i.e. where the Kraus operators of the 
corresponding map act only non-trivially on a finite, lo- 
calized number of subsystems. We have also considered 
lifetime of entanglement between subsystems of different 
size, which allowed us to determine the scaling behavior 
of entanglement under re-scaling of the size of the sub- 
systems. While for cluster states there is essentially no 
change in the scaling behavior with N, for GHZ states 
we found that (i) the lifetime of block-wise entanglement 
for any number of blocks that contain only a finite num- 
ber of particles to tends to zero as TV — > oo, while it can 
become finite if the blocks itself become macroscopic, i.e. 
TO ^ oo as ^ oo; (ii) the number of blocks M that 
remain entangled after a certain time t is finite, indepen- 
dent of the block size m. In addition, we have shown that 
for encoded entangled states the number of blocks that 
can be entangled after a certain time can be drastically 
increased and the lifetime of M-party entanglement can 
be enhanced. 

Our results suggest a remarkable robustness of certain 
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kinds of macroscopic entangled states -namely all graph 
states with constant degree- under various kinds of de- 
coherence. 
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Appendix A 

In this appendix we will give a more detailed analysis of 
the statement (ii) of Sec IV Dl Given an arbitrary channel 



(124) 



ij=0,l,2,3 



we first consider the question, whether it is possible to 
decompose it into a dcphasing channel T^p — ^^2' P + 
^—^a^pcTz with pz G [0, 1] followed by some arbitrary 
noise channel E' p = I]ij=o,i.2,3 9y fipcTj, i.e. 



£ = £' oV 



(125) 



Since the upper bound derived in Sec IV Dl becomes 
tighter with an decreasing dephasing parameter pz, we 
will also try to minimize Pz for those channels £, for 
which an extraction of a dephasing part is non trivial, 
i.e. Pz < 1- First of all, with the matrices P — (pij), 
Q (qtj) and 



M 




(126) 



we can rewrite Eq. H125|) as a matrix equation 



(127) 



which is linear in Q. For < < 1 it has the unique 
solution 



Q = El±lp + • P • M 



2pz 



2pz 



(128) 



kk' 



Since the matrices Q and P coincide with the states 
and cr|'^ obtained via the Jamiolkowski isomorphism (2^ 
(see Sec IV C|) . the conditions that Q actually corresponds 
to a completely positive and trace preserving map £' , are 
(a) trfc' Q = (Tq and (b) Q > 0, i.e. that Q is a density 
matrix. It is straightforward to show that tr^/ Q is the 
maximally mixed state ctq = ^1 on particle k whenever 
trfe/ P — (Tq and hence that for < < 1 condition 
(a) is always fulfilled. On the other hand the positiv- 
ity (b) imposes further constraints on the allowed pa- 
rameter range of pz, for which the channel £ permits 
a decomposition of the form (|125|l . For Pauli channels 



£p = J2i=o 1 2 3 Pi ^iP'^i the solution Q is again a diago- 
nal matrix: 

Q = T^-diag ( {pz + l)pQ + {pz - 1)P3, 



{Pz + l)Pl + {Pz - 1)P2, 

(Pz + l)P2 + {Pz - 

(Pz + 1)P3 + {Pz - l)Po ] 



(129) 



Now the positivity (b) yields the four inequalities: 

1 - ^ Po ^ l+Pz 



l+Pz 
l-Pz 

l+Pz 



P3 

< ^ < 
P2 



I -Pz 
l+Pz 

l-Pz 



(130) 



which imply that in particular either both pi and p2 [po 
and P3] are zero or none of them. Therefore in the case of 
a bitflip channel, as to be expected, no dephasing com- 
ponent can be extracted in the above sense. Moreover it 
is straightforward to see that the minimal value for pz, 
that allows for such an extraction is given by 



l-Pz 



< Qn 



Po Pi P3 P2 



l+Pz IP3 P2 Pa Pi 

1-9. 



(131) 



i-e. Pz ^ i+q^ ■ the trivial case of a dephasing chan- 
nel itself with parameter p, i.e. po = ^~2^^ Ps — 
and pi — P2 — 0, this clearly gives a minimal value 



_ l+3p 

4 ' 

we obtain the minimal value for 



for Pz — p. For the depolarizing channel (po — — 4 



Pi = P2 = P3 = ^ 
Pz = T+p^ ^^^^ ™ the case of the quantum optical chan- 
nel with /X = (po = 1 - (2p + g) , p 1 = p2 = P < i 
and P3 = q < \ ) one arrives at = 1 — — 



2p- 



Appendix B 

In this appendix we will give the proof to the state- 
ments (i), {a) and [in) of Sec IV El 

Proof of (i): Let C/i, U2, U3 denote the disjoint subsets 
on which a ax,cry or error occurs. From the stabilizer 
equations Kjf \G) = \G) it follows that a^\G) = crf^lG) 
and similarly (t^|G) = cr^<|G) = (Tf'=+'^|G). With this 
one obtains 



P = 



„|V|-|C/i|-|(72|-|(73l JUl\ \U2\ JC/3I 

/ . Po Pi P2 P3 



Ui,U2.U3CV 

'^f|G)(G|af , 



(132) 



where U = T{Ui + U2) + U2 + Us- With Qi = ^ the 
eigenvalues Xij therefore can be written as 



At/ = Po 



{Ui,U2,U3)eM{u) 



where 

Miu) = { (c/i, c/2, U3)\u,cv,u,n u, = 0, 
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u = r(Ui + f/2) + f/2 + [/3 } • 

This set may alternatively be written as 

Miu) = { (c/' \ [/" , u' n [/" , c/" \ c/') 
I [/',[/" c y , J7 = ru' + u" } 
= {{u'\ [Tu' + u),u'n {ru' + u) , 

{TU' + U)\U') \U' CV } . 

This is exactly the index set in Eq. JSUJ). 

Proof of (ii): We first note that for partial transposi- 
tion -^-^ generally 

= (Cl^)^ ® (D^^f- (C^n-^ ® C-f (133) 
holds. With this rule we can compute 

K^{\G){G\f^K^ = {K^\G){G\K^f^ = (|G>(G|f- , 

i.e. (|G)(G|)'^^ commutes with K'^ for aU fc G There- 
fore (|G)(G'|)^'' is again diagonal in the graph state basis 

{\U)g\U^ i-c. (|G)(G|)^- = Eacv A'^|C/)g(C/| ■ 
In order to determine the spectrum let us decompose U 
into (C/a, U A"^) accordingto the partitioning and use the 
Schmidt decomposition [g 



2^ 



with = (A',r^A'), the standard z-basis \G)^ = 

cr^jO)'^'^ on partition A and the graph basis \G)^'' — 
|G)q^^ on partition with respect to the pure graph 
state corresponding to the graph G\ A, which is obtained 
from G by removing all vertices in A. Now we can com- 
pute 



a;. := 



G{U\{\G){G\f^\U)G 



fA{Ai) + {UA,Al+A2) 



1=1,2,3,4 



{A,\A3){A2\A,)(T'A,\T'A2 + Uao) X 

x{r'As + UA^\r'A^) 
i\A\ 2^ ^ 



Ai ,A2QA: 
r'(Ai+A2) = Uj^a 

2\A\ >^ 



(134) 



Whereas Xjj vanishes for U CV with Ua'= ^ IniF', in the 
opposite case Ua" G ImF' we can choose an arbitrary Aq 
with T'Aq [/^c to simphfy (|134(l : 



81^1 



(-1) 



{UA:Ao) 



IkerT'l if C/^ € (kerT')"^, 
if f/A ^ (kerT)-^ 



Note that this is independent of the choice Aq, since a 
different choice Aq will differ from Aq only by an element 
A' e keiT', for which (Ua,A') = if C/a G (kerT')-^. 
Therefore one obtains the partial transpose for the pure 
graph state 

Ta Ikerr'l 



(|G)(G|)^- = 



21-41 



(-i)<^^^^>|x + r)G(x + r| (135) 



(kor r')^ X (Im r') 



The corresponding formula 1)83(1 for a general graph 
diagonal state can be deduced from Eq. ((135(1 by again 
using 1133 . 

Proof of (Hi): For Xu+k the estimation can be derived 
from Eq. I(8U(I : By adding [deleting] an element k to 
[from] the set TU' + U the corresponding sizes of the 
sets Ui = U'\ (TU' + U), U2 = U' n {TU' + U) and 
[/g = [TU' + U)\U' in the exponents of ^ can at 
most increase or decrease by one. Moreover, since Ui, U2 
and U3 are disjoint, the operation U t-^ U + k of adding 
or deleting k cannot simultaneously increase or decrease 
any two sets Ui. For example, if fc G Ui then 

\u'\{ru' + u + k)\ = [c/il-i, 
\u' n {TU' + u + k)\ = IC/2I + 1, 
\{ru' + u)\u'\ = IC/3I . 

In any case every addend 

\u'-{ru'+u+k)\ \u'n{ru'+u+k)\ \{ru'+u+k)-u'\ 
9i Q2 I3 

in l(8UII can be bounded from below by 5 x q^^^^ q]^^^ q\^^^ 
and from above by ^ x q!^'^^' q]^^^ q\^^^ since q < 1- This 
gives qXu < Xu+k < ^Xu- For Xu+n^ and Xu+N^+k a 
similar argument holds, if one rewrites 

U'CV 

\(u'+k)n{ru'+u)\ \(ru'+u)\(u'+k)\ 
X 92 la 



and 



Xu+Nt, + k - Po I1 



|((7'+/c)\(r;7'+;7-hfe)| 



U'CV 

\{u'+k)n(ru'+u+k)\ \(ru'+u+k)\{u'+k)\ 
X 1I2 Is 



In this representation we 'absorbed' Nk into the sum- 
mation index U' in both cases using TU' + U + Nk — 
T'(U' + k) + U. This concludes the proof of (ui). 
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